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Abstract. 

We study a two-dimensional array of coupled one-dimensional (ID) tubes of 
interacting bosons. Such systems can be produced by loading ultra-cold atoms in 
anisotropic optical lattices. We investigate the effects of coupling the tubes via hopping 
of the bosons (i.e. Josephson coupling). In the absence of a periodic potential along the 
tubes, or when such potential is incommensurate with the boson density, the system 
undergoes a transition from an array of incoherent Tomonaga-Luttinger liquids at high 
temperature to an anisotropic Bose-Einstein condensate (BEC), at low temperature. 
We determine the transition temperature and long wave-length excitations of the BEC. 
In addition to the usual gapless (Goldstone) mode found in standard superfluids, we 
also find a gapped mode associated with fluctuations of the amplitude of the order 
parameter. When a commensurate periodic potential is applied along the tubes, they 
can become ID Mott insulators. Intertube hopping leads to a dcconfincmcnt quantum 
phase transition between the ID Mott insulators and the anisotropic BEC. We also 
take into account the finite size of the gas tubes as realized in actual experiments. We 
map out the phase diagram of the quasi-lD lattice and compare our results with the 
existing experiments on such systems. 
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1. Introduction 

One dimension (ID) provides fertile grounds for studying the physics of strongly 
correlated quantum many-body systems. It is a well established theoretical result that, 
because of the enhanced role of quantum fluctuations in low dimensionality, there is no 
broken continuous symmetry and therefore long range order, even at zero temperature, 
in a ID quantum many body system. Instead, in a large class of systems possessing a 
gapless spectrum, power laws characterize the various correlation functions. This class 
of systems is known as a Tomonaga-Luttinger liquid (TLL). 

One of the simplest examples of a TLL is the Lieb-Liniger model [1 , which describes 
a system of non-relativistic bosons interacting via a contact (Dirac delta-function) 
potential. This model has been solved exaclty using the Bethe ansatz JT1. which has 
provided valuable insights into its thermodynamic properties and excitation spectrum. 
However, calculation of the asymptotics of correlation functions is not an easy task using 
this method, and in this regard it is nicely complemented by the harmonic-fluid approach 
introduced by Haldane pL The combination of these two methods has allowed us to 
understand that, as the the interaction is increased, the system exhibits no quantum 
phase transition. Instead, it smoothly crosses over from a weakly interacting (often 
called quasi-condensate [3 HE]) regime, where phase fluctuations decay very slowly, to 
a strongly interacting regime, the Tonks regime In the latter, repulsion between the 
bosons is so strong that it leads to an effective exclusion principle in position space, and 
makes the system resemble a gas of non-interacting spinless fermions in many properties. 
Thus it is sometimes stated that the bosons in the Tonks regime "fermionize" . 

This very simple model solved by Lieb and Liniger in 1963 has recently received a 
great deal of attention thanks to spectacular advances in confining quantum degenerate 
gases of alkali atoms in low dimensional traps [7J. While at the beginning only weakly 
interacting ID Bose gases were available jSHH], the application of a deep optical potential 
parallel to the ID axis to arrays of ID Bose gases has recently allowed the group at 
ETH (Zurich) ^21 IH| to demonstrate a quantum phase transition from a superfluid 
TLL to a ID Mott insulator, which is the low dimensional analog of the superfluid to 
Mott insulator transition first demonstrated in the 3D optical lattice by Greiner and 
coworkers ^2] ■ Furthermore, Kinoshita et al. (THj have succeeded in reaching the Tonks 
regime by creating an array of tightly confined ID Bose gases. A different path ^] 
to realize a strongly interacting ID Bose gases has been followed by Paredes et al. [16J, 
who applied a deep periodic potential to an array of (otherwise weakly interacting) ID 
Bose gases. Thus they were able to dramatically increase the ratio of the interaction to 
the kinetic energy of the atoms. In both cases [HUGE!, the observation of fermionization 
signatures was reported. 

In the past, much experimental effort focusing on low dimensional systems has 
been spent on studying electronic systems such as spin chains and quasi- ID metals (for 
a review see jlZj). Theoretically at least, electrons in ID become a TLL (provided there 
are no gap-opening perturbations) as soon as interactions are taken into account in a 
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non-perturbative way. However, for the electronic systems comparison between theory 
and experiment has proved problematic partly because of material-specific issues such as 
the existence of many competing phenomena at similar energy scales (e.g. in Bechgaard 
salts 0111) or the always-present disorder in solid state systems [T7]. Furthermore, 
one major complication stems from the fact that these materials are indeed three 
dimensional (3D), as there is always a small amplitude for electrons to hop from chain to 
chain, whereas Coulomb interactions (often poorly screened because of bad metallicity) 
also enhance 3D ordering tendencies. These inter-chain couplings cannot be neglected 
at low temperatures, which leads to the interesting physics of dimensional crossover: 
at temperatures high compared to the energy scale determined by these inter-chain 
couplings, the system effectively behaves as a collection of isolated ID chains. However, 
when the temperature is lowered below this scale, 3D coherence develops between the 
chains and, at zero temperature, a complete 3D description of the system is required. 
Much exciting physics is expected to occur in this crossover from the exotic ID TLL 
state to a conventional 3D (but strongly anisotropic) Fermi liquid, or an ordered 3D 
phase [T7j. Indeed, these phenomena may well have relevance to the pseudogap (and 
other anomalous) behavior of the quasi-2D high temperature superconductor materials. 
For instance, one of the many interesting questions that arise in this context is the 
following: how [19J do the exotic low-energy excitations (bosonic collective modes) of 
the ID system transform into the Landau quasi-particles characteristic of Fermi liquids 
and which resemble more the constituent electrons? 

The dimensional crossover phenomena described above has also a counterpart in 
Bose gases loaded in quasi-lD optical lattices. From a theoretical point of view, the fact 
that bosons can undergo Bose-Einstein condensation and, therefore, can be collectively 
described using a macroscopic quantum field makes more amenable these systems to a 
theoretical analysis relying on a mean field approximation. Such a mean-field analysis 
is undertaken in this work: here we treat the intertube hopping via a mean field 
approximation, while quantum fluctuations within a tube are treated accurately at low 
energies using the harmonic-fluid approach. 

From the experimental point of view, the recent availability of gases of ultra- 
cold bosonic atoms in optical lattices has revived the interest in understanding the 
properties of coupled systems of bosons in ID. Furthermore, the already demonstrated 
high degree of tunability of these systems can provide a clean setup to study the physics 
of dimensional crossover. Thus, as we show in sectionElof this work, the phenomenology 
of Bose-Einstein condensation in the most anisotropic version of these lattices is very 
different from the one exhibited by Bose gases condensing in the absence of a lattice or 
even in shallow isotropic optical lattices. 

In addition to studying the properties of Bose condensed systems in lattices, some 
experimental groups have applied a tunable optical potential along the longitudinal 
direction of an array of effectively uncoupled ID Bose gases. It has been thus 
possible as mentioned above, to drive a quantum phase transition from the ID 
superfluid (TLL) state to a ID Mott insulator. Under these circumstances, the effect 
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of a small hopping amplitude (i.e. Josephson coupling) between neighboring tubes in 
the array becomes most interesting: as we show in section HJ for a finite value of the 
hopping amplitude, the system undergoes a new type of quantum phase transition, 
known as "deconfinement" transition, from the ID Mott insulator to an anisotropic 
Bose-Einstein condensate. Thus the tendency of bosons to develop phase coherence 
overcomes the localization effect of the optical potential. As explained in section El we 
believe that some features of this deconfinement transition may have been observed in 
the experiments of Ref. 

This paper is organized as follows: section El defines our model of a 2D array of ID 
tubes of bosons coupled by weak intertube tunneling, and describes the harmonic-fluid 
approach that allows us to derive a low-energy effective Hamiltonian. Then we discuss 
separately the phenomenon of Bose-Einstein condensation and the properties of the 
condensate in an anisotropic lattice (section EJ), and the deconfinement quantum phase 
transition (section In the former case, we obtain the zero-temperature condensate 
fraction, the condensation temperature, and the excitations of the condensate using the 
mean field approximation and including gaussian fluctuations about the mean field state. 
We also show that some of the results obtained by these methods are recovered from the 
self-consistent harmonic approximation ( section l3.2|) . These calculations are described 
in pedagogical detail in Appendix |Appendix B[ and Appendix C| Moreover, in 
section IrS.l.rSI we discuss the excitation spectrum of the Bose-condensed phase, and in 
particular the existence of a gapped mode related to oscillations of the amplitude of 
the condensate fraction, in the light of a phenomenological Ginzburg-Landau theory. 
Differences between this theory and the more conventional Gross-Pitaevskii theory, 
where the gapped mode is absent, are exhibited. The transition from the (anisotropic) 
regime described by the Ginzburg-Landau theory to the one where Gross-Pitaeskii 
theory holds is also outlined there. 

The deconfinement transition is described in section 0] There we show how 
the renormalization group can be used to understand the competition between 
the localization in the optical potential and a small tunneling amplitude between 
neighboring ID Bose systems. To assess the nature of the deconfinement transition, 
we map (see section 14.2)1 the mean-field Hamiltonian, at a specific value of the system's 
parameters, to a spin-chain model. Finally, section discusses the effect of the finite size 
of the ID system that are realized in actual experiments, and how it modifies the phase 
diagram obtained from the renormalization-group analysis. We argue that, for an array 
of finite ID Bose gases, there is a extra energy scale below which the tunneling of atoms 
between tubes is effectively blocked. Finally, in section |H1 we present a discussion of the 
experimental observations and how they compare with the predictions of our theory. A 
brief account of some of our results has appeared elsewhere 
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2. Models and Bosonization 

In this section we introduce the models for an interacting Bose gas in a quasi-lD optical 
lattice that we study in the rest of the paper. We also briefly describe the harmonic- 
fluid approach (henceforth referred to as "bosonization"), which allows us to obtain a 
convenient effective low-temperature description of these models. 

2.1. Models of a Bose gas in a quasi-lD optical lattice 

In this work we study the low-temperature equilibrium properties of systems recently 
realized by several experimental groups [TJ EH CHI I2H 1221 I2H] which have succeeded 
in loading ultracold atomic gases in highly anisotropic (quasi-lD) optical lattices. In the 
experiments carried out by these groups, a Bose-Einstein condensate of ultracold alkali 
atoms (in particular, 87 Rb) is placed in a region of space where two pairs of non-coherent 
and mutually orthogonal counter-propagating lasers intersect. The intensity of the lasers 
is increased adiabatically, and because of the AC stark effect, the atoms experience 
a periodic potential |21j V±(y,z) = Vq v sin 2 {2ny/\) + Vq z sin 2 (2irz/\), where usually 
Voy — Voz — Voj_, with Vo_l oc Eq, E being the maximum intensity of the electric field of 
a laser whose wavelength is A. The strength of the optical potential, Vq±, is frequently 
given in units of the recoil energy, Er = h 2 7r 2 /(2Ma 2 ), where a = A/2 is the period of 
V±(y, z ) an d M the atom mass f. For V ± fi-^D, where ~ 47r/i 2 a s p iD / M j2S] is 
the chemical potential of the BEC (a s being the s-wave scattering length of the atom- 
atom interaction potential, and pjj D the mean atom density), the gas atoms are mainly 
confined to the minima of V±(x, y), thus forming long, effectively ID, gas tubes as shown 
in figure [*[] (there are additional harmonic potentials along the x, y, z-axes that trap the 
atoms). Under these conditions, we can project the boson field operator ty(x, y, z) onto 
the lowest (Bloch) band of the periodic potential, 

*(x,y,z)^Y, W *-(y> z ) iS n( x )> W 

R 

where R = (m y , m z ) a, with integers = —Mi/2, . . . , +Mj/2 (i — y,z), is the location 
of the different potential minima; Wn(x,y) is the Wannier function of the lowest Bloch 
band centered around R. The system just described is called a 2D (square) optical 
lattice. In typical experiments M y ~ M z ~ 10 2 , and thus the number of gas tubes is of 
the order of a few thousand separated half a laser of wave-length, A/2 ~ 0.4 /im. 

In addition, in several of the experiments [TU1 1221 EI] an extra pair of lasers 
was applied to define a longitudinal periodic potential V\\(x) = Vq\\ sin 2 (2tix/\) = 
| Vo| | +u(x), and u(x) = u cos(Gx), where uq = — Vo||/2 and G = An/X = 2n/a. Since, 
for a commensurate potential, z.e.when the period a matches the mean interparticle 
distance, the variation of the strength of u(x) is responsible for the transition [2*1 IT7] 
to a Mott insulating state where the bosons are localized, we shall henceforth refer to 

| For 87 Rb, Er/% = a/A 2 , where a ~ 14.424 ^m 2 kHz and A is the laser wavelength in /im. Thus, for 
A = 0.826 /mi used in the experiments of Ref. !() . Er/H = 21.14kHz 
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Figure 1. Sketch of a two-dimensional array of coupled tubes (also called quasi- ID 
optical lattice in this work). Bosons can hop from tube to tube with amplitude J, 
which can be controlled by the height of an optical potential transverse to the tube 
axis. Along the tube, a periodic optical potential (for clarity, shown here for only one 
of the tubes) can be also applied. This potential can drive a quantum phase transition 
where each tube transforms into a ID Mott insulator. 



u(x) as the "Mott potential". In this work we are interested in the extreme anisotropic 
case where the strength of this potential |wo| ~ Vo\\ <C Vo±. In experiments, this was 
achieved by increasing Vq± so that the motion of the atoms effectively becomes confined 
to a one-dimensional tube, where they undergo only zero-point motion in the transverse 
directions. This picture is equivalent to the projection onto the lowest Bloch band 
described above. Just like in three dimensional space, the atomic interactions in the 
tube are still well described by a Dirac-delta pseudo-potential, v a (x) = g\o 8(x), where 
the effective coupling has been obtained in Ref. [27J: 



9m 



(2) 



(1 -Ca s /V2£ X )M£ 2 ± ' 

The constant C ~ 1.4603, and the oscillator length in the transverse direction, £j_ = 
(/i/(Mwi))^ 2 (where huj±/E R ~ 2(V ±/ 'E R ) l l 2 for deep lattices [TJ], is the transverse 
oscillation frequency). 
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Taking into account the above considerations, we are led to consider the following 
Hamiltonian: 

H = 22 dx WV7 d ^k( x ) d x^n(x) + u(x)pr(x) + - / dx' v a (x - x')p n (x)pn(x') 
^ Jo \_2M 2 Jo 

~i E [ L dx[^ R (x)^ Il/ {x)+h.c.] , (3) 

A (R,R'> J0 

where Pr(x) = ^^(z^r^) is the density operator at point x and lattice site R; the 
field operator \I/r(x) obeys [\I/r(x), ^^/(V)] = ^rr' 5(x — x'), and commutes otherwise 
as corresponds to bosons. The last term in equation (jHJ) describes the hopping of bosons 
between two neighboring tubes. This term is also known as Josephson coupling. The 
hopping amplitude 

e _ 2(W ^ ; (4) 

for Vo_l ^> -Ejj. We shall assume a very deep transverse potential so that J where 
/Wid = ^id(Po) is the chemical potential in each tube, as can be obtained from the exact 
solution of Lieb and Liniger |X| (or measured experimentally). Since experimentally 
Hid < 0.1 Tioj± |7 3 , assuming J pm requires at least that Vo± > 10 Er. Such limit 
allows us to regard the system as an array of weakly coupled ID interacting Bose gases. 

In the experimental systems, the Mott potential u(x), besides the periodic part 
described above, also contains a slowly varying harmonic term that confines the atoms 
longitudinally. Additionally, there are also harmonic traps in the transverse directions. 
This makes the array, not only finite, but also inhomogeneous. Whereas we shall discuss 
finite-size effects further below (see section |SJ), it has proven so far intractable in ID 
(except for a few limiting cases to deal with the harmonic confinement analytically 
and only numerical results are available [SHI EU EH E3 EH EH] • Thus, with the exception 
of sect ion [7| where we discuss the experimental consequences of our work, we shall neglect 
the effect of the harmonic potential, and treat the system as an array of gas tubes of 
length L and hard- wall box (open) boundary conditions [HSIEZIEH], each one containing 
iVo = PqL bosons (in typical experiments, iVo ranges from a few tens to a few thousand 
particles). This makes sense, however, as much of the interest of this work is focused on 
computing the phase diagram as well as some of the thermodynamic properties of the 
different phases, for which one is required to take the thermodynamic limit L — > oo at 
constant density p = N /L. 

Finally, when the axial potential u(x) is made so deep that |u | ^> pm (but still 
l M o| <C Vo_i_), it becomes convenient to perform a further projection of ^r^) onto the 
lowest Bloch band of the Mott potential u{x). Thus, 

M 

*r0*0 ^ w m{x) b m n, (5) 

m=l 
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where [6 m R, b m i R ] = Sr.^ 8 mtm i , commuting otherwise, and w m (x) is the Wannier orbital 
centered around site m. This projection leads to the so-called Bose-Hubbard model [39J: 

H BU =J2Hn + Hj, (6) 

R 

#R = -TT E & LR & »«'R+7E n mR» ( 7 ) 

Z (m,m'> ™ 

= " y E Kr^r' + h.c] , (8) 

(R,R'),m 

where n m R = ^R^mR is the boson number operator at site (m, R). Note that \u \ <C V ± 
implies that J x ^> J, and therefore we are dealing with a very anisotropic (quasi-lD) 
version of the Bose-Hubbard model. 



2.2. Bosonization and low-energy effective theory 

In order to treat the interactions in a non-perturbative way, it is convenient to rewrite (|HJ) 
in terms of collective variables that describe phase and density fluctuations. The method 
of bosonization H3 HOI E±| allows us to obtain such a description. The description is 
accurate as long as one is interested in the low-temperature and low-frequency properties 
of the system. In this subsection we shall briefly outline this method, but we refer the 
reader to the vast literature on the subject HZ1 EH EH] for more technical details. 

To study the model of equation Q, we use the density-phase representation of the 
field operator. To leading order, 

*r(*) = (AbPo) 1/2 exp [i9 R (x)) , (9) 

where is Ab is a model dependent (i.e. non-universal) prefactor (see below). The density 
operator, 



Po d x (f) K (x) 



71 



e 2im{iTp x-(l> R (x)} (10) 



The fields 6*r(x) and 0r(x) describe collective fluctuations (phonons) of the phase and 
the density, respectively §. Since the density and the phase are canonically conjugate 
variables, it is possible to show |H3 Ej that ^r(x) and 0r(x) obey the commutation 
relation: [0 R (x), d x i6-Ri(x')] = «5rr/5(x — x'). Furthermore, it is assumed that <9 x 0r(x) 
is small compared to the equilibrium density po- Note that this approach does not 
assume the existence of a condensate and thus properly takes into account the quantum 
and thermal fluctuations, which are dominant in ID. Indeed, it is the particle density 
and not a "condensate" density that appears in the density-phase representation of 
equations ()9ll0p (unlike usual mean field approaches). Furthermore, the terms of 
equation with m/fl account for the discrete particle nature of the bosons, which 
is crucial for obtaining a correct description for the transition from the ID superfluid 
phase to the ID Mott insulator [2| 126) ITT] HTj EE2] . 

§ We follow here the notations of ^Ij ■ Note that compared to [371 |23j 4> an d & are exchanged and the 
density field, <fi, in this work differs by a minus sign from the density field, 9, in those references. 
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Using equations JHJ) and (|1U|). the Hamiltonian in (jSJ) becomes 

hv "~ rL ' ^K{d x d^{x)f + ^{d x Mx)f 



2n ^ Jo 



ZilT CL J 

rL 

V / cli cos(fe(i)-^(i)), (11] 



R 

2 ™ 2 ( R . R ' 



where 5 = (G — 2ttpq)/tt = 2(a _1 — po) is the mismatch between the density po an d the 
periodicity of the Mott potential. The two last terms are related to the Mott potential 
and the Josephson coupling, respectively. We use the following dimensionless couplings 
to characterize their strength: 

27rJ(p a ) 2 2iru (p a ) 2 

9J = — r A B , 9u = z (12) 

hv s po nv s po 

where a ~ hv s /p,iD is a short-distance cut-off. 

The dimensionless parameter K and the sound velocity v s depend both on the 

density, po, and on the strength of the boson-boson interaction. Indeed, for the Lieb- 

Liniger model, they depend on a single dimensionless parameter, the gas parameter pQ 

7 = Mgi£)/h 2 po- Analytical results for these parameters are only available in the small 

and large 7 limit (see ref. jHZ] for results in the intermediate regime): 

VsK" 1 ~ v F [l - 87- 1 ] (7>1), (13) 

v s K' 1 ~ w F 7/vr 2 (7 < 1), (14) 

whereas Galilean invariance fixes the product |2"1 137|: 

v hn P° ft z\ 

V S K = V F = —-r- 15 

M y ' 

Therefore, 1 < K < +00 for the Lieb-Liniger model, K = +00 corresponding to free 
bosons and K = 1 to the Tonks limit. For the Lieb-Liniger model, the formula [SZj: 

As{K)^(-) (16) 



7T / 

provides a reasonably good interpolation (accurate to within 10% of the exact results) 
for the non-universal prefactor of the field operator. 

It is worth stressing that in (fTTj) we have retained only terms which can become 
dominant at low temperatures. In the language of the rernormalization group (RG) 
these correspond to marginal or relevant operators. Whenever, for reasons to be given 
below, any of the terms retained in equation (jllj) becomes irrelevant in the RG sense, 
we shall drop it and consider the remaining terms only. We emphasize that the validity 
of (fTTj) is restricted to energy scales smaller than min{T, pim}, and therefore its analysis 
will allow us to determine the nature of the ground state and low-energy excitations of 
the system. 

Interestingly, the effective Hamiltonian of equation (fTTj) also describes the low 
energy properties of the anisotropic Bose-Hubbard model of (JBJ)- This result can be 
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derived by a lattice version of the bosonization procedure ^7j, followed by to 
the continuum limit. However, the relationship of the effective parameters K, v s , g u , 
and gj to microscopic parameters J x , J, and U is no longer easy to obtain analytically 
(except for U small and certain limits in ID ^3]). Thus, in this case K, v s , g u , and 
gj must be regarded as phenomenological parameters, which must be extracted either 
from the experiment or from a numerical calculation. 

3. Bose-Einstein Condensation 

Let us first discuss the phenomenon of Bose-Einstein condensation in anisotropic (quasi- 
1D) optical lattices. This is the phase transition from a normal Bose gas to a Bose- 
Einstein condensate (BEC), which exhibits some degree of spatial coherence throughout 
the lattice. Coherence arises thanks to the Josephson term, which couples the ID 
interacting Bose gases together. In the absence of the Josephson coupling and the 
Mott potential (i.e. for J = and u(x) = 0), the system behaves as an array of 
independent TLL's 011311011211: the excitations are ID sound waves (phonons) and, 
at zero temperature and for L — ► oo, all correlations decay asymptotically as power- 
laws (e.g. phase correlations (e ieR ^e~ l9li '^) ~ <5r,r' |po^| _w f° r x ^ Po 1 ), which 
implies the absence of long-range order. However, as soon as the Josephson coupling 
is turned on, bosons can gain kinetic energy in the transverse directions, and therefore 
become delocalized in more than one tube. This process helps to build phase coherence 
throughout the lattice and leads to the formation of a BEC. Nevertheless, the building of 
coherence just described can be suppressed by two kinds of phenomena: one is quantum 
and thermal fluctuations, which prevent the bosons that hop between different gas 
tubes from remaining phase coherent. The other is the existence of a Mott potential 
that is commensurate with the boson density. For sufficient repulsion between the 
bosons, this leads to the stabilization of a ID Mott insulator, where bosons become 
localized. Whereas we shall deal with the latter phenomenon in section we deal with 
the fluctuations in this section. 

In this section we have in mind two different experimental situations: one is a 
2D optical lattice where the Mott potential is (experimentally) absent. The other is 
when the Mott potential is irrelevant in the RG sense. The precise conditions for this 
to hold will be given in section HI but, by inspection of JUJ), we can anticipate that 
one particular such case is when the Mott potential is incommensurate with the boson 
density (i.e. 5 ^ 0). Under these circumstances, the cosine term proportional to g u 
oscillates rapidly in space and therefore its effect on the ground state and the long 
wave-length excitations is averaged out. 

3.1. Mean field theory 

To reduce the Hamiltonian (jSJ) to a tractable problem, we treat the Josephson coupling 
in a mean-field (MF) approximation where the boson field operator ^r(o;) is replaced 
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by ip c + 5\E'r(x), where ipc = (^r( x )) an d 8^r(%) = ^n(%) — 4>c- We then neglect the 
term oc 5^h(x)5^/-r>(x), which describes the effect of fluctuations about the mean field 
state characterized by the order parameter ip c . This approximation decouples the ID 
systems at lattice sites R 7^ R' so that the problem reduces to a single "effective" ID 
system in the presence of a (Weiss) field proportional to \ip c \- The bosonized form of 
the Hamiltonian for this system is (we drop the lattice index R): 

K (dj(x)) 2 + K- 1 (d x( p(x)f 

- 2Jz C \j 'AbPo\?Pc\ f o dx cos0(x) + Jz c L\?p c \ 2 , (17) 

where zq is the coordination number of the lattice (zc = 4 for the square lattices 
realized thus far in the experiments). The phase of the order parameter 9 C = arg-0 c 
has been eliminated by performing a (global) gauge transformation 8(x) — > 6{x) + 6 C , 
and henceforth we shall take ip c = \ip c \ to be real. We note that, up to the last 
constant term, H^ F defines an exactly solvable model, the sine-Gordon model, of which 
a great deal is known from the theory of integrable systems as well as from various 
approximation schemes |17H IT?) HH] . The mean-field Hamiltonian must be supplemented 
by a self-consistency condition, which relates the effective ID system to the original 
lattice problem: 

^ c = (*n(aO) = 2 (A BPo ) 1/2 (cos#(x))mf (18) 
Note that the average over the cosine must be computed using the eigenstates of H^ F . 



H, 



MF 
eff 
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3.1.1. Condensate fraction at T = 0: To obtain the condensate fraction at zero 
temperature, we first compute the ground-state energy of H^ F using the following result 
for the ground state energy density of the sine-Gordon (sG) model 



BsG{\^c\) 



tan 



Tip 



(19) 



4hv s " " 2 ' 

where p = (3 2 / (8n - (3 2 ) = l/(8K-l), with f3 2 = n/K. The excitations of the sG model 
are gapped solitons and anti-solitons (and bound states of them for (3 2 < An) [4T?1 ITTj. 
The soliton energy gap (also called soliton "mass") is [43J: 

-, p+i 

m s (\ip c \) 



A,(|^c|) = hv. 



and 



m a (\if>e\) 



«(?) 



k(p) 
zcJ(AbPo) 1 / 2 

hv s OQ 



2(p+l)/p 



\4>c 



ir 



p+i 



7T Y 



1 

p+1 



2+1 
2 



2r(f) 



(20) 



(21) 



(22) 



where T(x) is the gamma function. Hence, the ground state energy of is given by: 
E™ F (\ij c \)/L = £ sG (\ij c \) + Jz c \^ c \ 2 . (23) 
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The self-consistent condition ()18|) is equivalent to minimizing the mean-field free energy 
with respect to the order parameter. At T = 0, the entropy is zero and the condensate 
fraction is thus obtained by minimizing E^ F (|^ C |) with respect to ip c = \ip c \, which 
yields: 



In the above expression, we have introduced the dimensionless parameter r] = 
AB{K){pQa ) l / 2K . We note that the condensate fraction at T = has a power-law 
dependence on the amplitude of the Josephson coupling, i.e. 



The dependence of the exponent on the parameter K, which varies continuously with 
the strength of the interactions between the bosons in the tubes, is a consequence of 
the fact that the appearance of a BEC at T = 0, as soon as J is made different from 
zero, is affected by the characteristic ID quantum fluctuations of the TLL's. This 
behavior is markedly different from what can be obtained by a Gross-Pitaevskii mean- 
field treatment of the system that overlooks the importance of ID quantum fluctuations. 
Further below we shall also encounter other instances where the present results strongly 
deviate from Gross-Pitaevskii theory. 

To show the effectiveness of the fluctuations to cause depletion of the BEC, let 
us estimate the fraction of bosons in the condensate at T = using equation 
For K — 1 (i.e. the Tonks limit, where phase fluctuations are the strongest), we take 
a p = 1 and assume that pm = 0.1 hu± [7j. Therefore, for V ± = 20 E R , the hopping 
amplitude J ~ 3 x 10~ 3 fi 1D , and hence, using equation (J23|) . Po Vc ~ 20 %. We also 
note that for weakly interacting bosons, i.e. as K — > +oo, ?/r? — > p , as expected. 

In the experiments of Ref. ^U], very small coherence fractions were observed. The 
experimental coherence fraction is defined as the fraction of the total number of atoms 
that contributes to the peak around k = in the momentum distribution, as measured 
by time of flight. The fraction was quantified by fitting a gaussian distribution to the 
peak ^U|. We note that such a procedure yields a non- vanishing coherence fraction even 
for the case of uncoupled ID tubes, where the momentum distribution exhibits a peak 
with a power-law tail ~ \k\ 2K~ l and a width ~ max{L, hv s /T} jSZHSEj- Thus, identifying 
the condensate fraction with the experimental coherence fraction may not be appropriate 
as there is no condensate in such a ID limit. Furthermore, as we show below, in the 
anisotropic BEC phase discussed above, this experimental procedure would also pick up 
a non-condensate contribution to the coherence fraction. Nevertheless, provided heating 
effects and thermal depletion of the condensate are not too strong, one can regard the 
experimental coherence fraction as an upper bound to the condensate fraction, and 
since a strong decrease of the coherence fraction was reported in the experiments of 
Ref. JU] j our results showing a strong depletion of the condensate are consistent with 
the experimental findings. 




(24) 




(25) 
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3.1.2. Condensation temperature: At finite temperatures, thermal (besides quantum) 
fluctuations cause depletion of the condensate. Ultimately, at a critical temperature 
T = T c , the condensate is destroyed by the fluctuations. To compute the critical 
temperature we follow Efetov and Larkin JTIj and, taking into account that near T c the 
order parameter ip c {T) is small, consider a perturbative expansion of the self-consistency 
condition (fTSj) in powers of ip c . To lowest order, 



^ c (T) = (*(x,t)> 



(T{*(a;,T)e-H^'*w(rOj >o 



mf — nw 



n/3 

JzcMT) I dx' I dT' gi (x',T')+0(tf). (26) 



where \l/(x, r) = (poA B ) 1/2 e ie ^ T \ t being the imaginary time, T the ordering symbol 
in (imaginary) time, and the Weiss field, 

H w = -JzcMT) [ L dx [tf (x) + h.c] . (27) 
Jo 

The correlation function: 

9i(x, r) = -~(T [V(x, r)¥(0, 0)]) , (28) 

where the average (. . .) is performed using the quadratic part of H^ F . Hence, 
analytically continuing to real time, equation (}2l)j) reduces to: 

1 + Jz c g?(q = 0,u = 0;T c ) = 0, (29) 

which determines T c ; g^(q,u>) is Fourier transform of the retarded version of (|28|). This 
function is computed in the |Appendix B| (cf. equation ()B.15|) ). and leads to the result: 

\hv s po) FK \hv s po)' 



( \ f 1 - i 



where 

F(A-)^M B (A-) 8i „ ( -j B ^-,l--j, (31) 

with B(x, y) = T(x)T(y)/T(x, y) the beta function. Thus, the dependence of T c on J is a 
power-law whose exponent is determined by K. It is interesting to analyze the behavior 
of T c in the limit of weak interactions, i. e. for K S> 1. Using that F(K ^> 1) ~ K, along 
with Galilean invariance, equation (fTKjl . v s = vf/K, we find that 

X\ 2 (J 



r w (32) 

where e F = h7ip vp/2 = (hnp ) 2 /2M. Thus, we see that as K — > +oo, T c ~ J 1 / 2 . This 
result comes with two caveats. The first one is that, as we turn off the interactions 
(i.e. for giD — ► 0) so that K can diverge, the chemical potential pm — gmPo —> 0. 
Therefore since the effective Hamiltonian, equation (jllj) . is valid only for energy scales 
below p 1D , the range of validity of the result is going to lower and lower J in this 
limit. Indeed one cannot simply let K — > oo and keep J fixed. One would get 
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T c / H\d ~ (cf / ' I^id) 1 ^ 2 ^ 'I Vid) 1//2 —* +00. On the other hand, we know that as gm — > 
T c must remain finite because a non-interacting Bose gas in a 2D optical lattice indeed 
undergoes Bose-Einstein condensation (see below). Therefore, if T c j ' \i\r> diverges as 
A 4 id — * for K ^> 1, it is because of this fact. The second caveat is that the behavior 
T c ~ J 1 / 2 implied by (jH2| is not the correct one in the non- interacting limit. To see 
this, consider the condition for the Bose-Einstein condensation of a non-interacting 
Bose gas in a 2D optical lattice (we take L — > 00 but keep the number of lattice sites 
M = M y x M z finite): 

Vol 

where e(q) and e_i_(Q) = —2 J (cos Q^a + cos Q z a — 2) are the longitudinal and transverse 
dispersion of the bosons, respectively. In the latter case, we have chosen the zero-point 
energy (i.e. the bottom of the band) to be zero. In the J —>■ limit the ID systems 
become decoupled and T c vanishes (i.e. there is no BEC in ID). Therefore, we expect 
T c to be small for small J. It is thus justified to expand about Q = (0, 0) the transverse 
dispersion e_i_(Q) — Ja 2 Q 2 . Furthermore, e(q) = h 2 q 2 /2M, for free bosons, and thus 
(|3*3*|) leads to the law T c ~ J 2 / 3 , which is different from (j32j) . The difference stems 
from the assumption of a linear longitudinal dispersion of the excitations implied by 
the quadratic part of (|TT|) and (JT7J). Had we insisted in assuming a linear longitudinal 
dispersion of the bosons e(q) = hv \q\, where v has dimensions of velocity, we would 
have found that T c ~ J 1 / 2 as implied by (|32)l . This argument implicitly assumes that in 
the weakly interacting limit, the difference between particles and excitations vanishes. 
This seems quite reasonable, as the spectrum of excitations of the non-interacting Bose 
gas coincides with the single-particle spectrum. 

Nevertheless, we can expect equation (J30|) to provide a good estimate of the 
condensation temperature as long as J/fim "C 1. This is of course particularly true 
in the strongly interacting limit where K ~ 1. The interpolation formula for the 
prefactor Ab, equation (jlfij) . allows to be more quantitative than Efetov and Larkin 
and to obtain an estimate of T c for K = 2, for instance. For the Lieb-Liniger model this 
value of K corresponds to 7 ~ 3.5 jHZ| and to /xidAf — 0.367. For J ~ 3 x 10 3 /iiD 
(i.e. Voj_ = 20 E R ), equation (jSUj) yields T c /a*id — 0.6 and therefore T c /h ~ 10 kHz or 
T c ~ 70 nK. We note this temperature seems to be well above the experimental estimate 
for the temperature of the cloud, T < 10~ 3 hu± w 10~ 2 [j, 1D Thus, we expect the 
experimentally realized systems to be in the anisotropic BEC phase described in this 
section, and the optical lattice to exhibit 3D coherence. 

3.1.3. Excitations of the BEC phase: Let us next compute the excitation spectrum 
of the BEC. Due to the existence of coherence between the ID systems of the array, 
this turns out to be very different from the "normal" Bose gas (TLL) phase, where 
the gas tubes are effectively decoupled (i.e.ip c = 0) because fluctuations destroy the 
long range order. The most direct way of finding the excitation spectrum is to consider 
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configurations where the order parameter varies slowly with x and R and time. Thus 
the order parameter becomes a quantum field, ^/ c (x.R, t). One then asks what is the 
effective Lagragian, or Ginzburg-Landau (GL), functional, £gl(^, R, t), that describes 
the dynamics of \l/ c (x, R, t). There are at least two ways of finding this functional: one 
is to perform a Hubbard-Stratonovich transformation and to integrate out the original 
bosons to find an effective action for the fluctuations of the order parameter about a 
uniform configuration. To gaussian order (random phase approximation, RPA) this is 
carried out in |Appendix A| This method has the advantage of allowing us to make 
contact with the microscopic theory (the Hamiltonian (jllj) in our case) and thus to 
relate the parameters that determine the dispersion of the excitation modes to those 
of the microscopic theory. Another method, which we shall follow in this section, is to 
guess £gl using symmetry principles only. To this end, we first write Cql = T — V, 
where the kinetic energy term T depends on gradients of \l/ c (x, R, t), whereas V depends 
only on \l/ c (x, R, x, t). Both T and V are subject to invariance under global Gauge 
transformations: \l/ c (x, R, t) — > e t9 ° \l/ c (x,R, t). 

Using global gauge invariance, to the lowest non-trivial order in \l/ c , we can write 
V = a|\l/ c (x, R, t)\ 2 + b\^/ c (x, R, t)| 4 . Since we are interested in the ordered phase where 
— V'o we must have a < 0. Thus, we can write V as follows: 

V = ^(|* c (x,R,t)| 2 -|^ c | 2 ) 2 , (34) 

where A is a phenomenological parameter. On the other hand, in order to guess the form 
of the kinetic term T global gauge invariance alone is not sufficient. Further insight can 
be obtained by close inspection of the effective low-energy theory described by (jllj) . If 
we write it down in Lagragian form: 



hK 



1 



hv s gj 



£ eS ( x ,t) = —Y, -(d t 9 K y -v s (d x e R Y + T ^ r £ cos(0 r -0rO, (35) 

Z7T R 



27rO (R,R'> 




(36) 



we see that it is invariant under Lorentz transformations of the coordinates (x, t) where: 

/ v s t' \ ( cosh (3 sinh j3 
\ x' J \ sinh (3 cosh (3 

being f3 a real number. Thus, T can only contain Lorentz invariant combinations 
of \d t ^ c (x, R, t)\ 2 , \d x $ c (x, R, t)\ 2 , and |Vr^ c (x, R, t)| 2 , which are the lowest order 
derivative terms allowed by global gauge invariance. This leads to 

T(x,R,t) = [\d t V c (x,R,t)\ 2 -v 2 \d x V c (x,R,t)\ 2 -v 2 ± \V^ c (x,R,t)\ 2 } , (37) 

where vn = v s and Z<i and v± are again phenomenological parameters. Higher order 
terms and derivatives are also allowed, but the above terms are the ones that dominate 
at long wave-lengths and low frequencies. Therefore, the GL functional reads: 

HZ 

C GL (x,R,t) = -^[\d t ^ c (x, R,t)\ 2 -v^\d x q c (x, R,t)\ 2 

-vl\V R y c (x,R,t)\ 2 } - — (|M/ c (x,R,t)| 2 -|^ c | 2 ) . (38) 
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We find the excitation spectrum by setting ty c (x, R, t) = [ip c + Vc( x , R<> t)} e iec( - x ' R '^ 
and keeping only terms up to quadratic order, which yields: 

Cgl(x, R, t) ~ ^ [(d t r] c ) 2 - v 2 (d x r] c ) 2 - v\ (V R r/ c ) 2 - Z^X^r] 2 

+ [{d t e c f - v\ (d x 6 c ) 2 - vl (V R # C ) 2 ] (39) 

Hence cf 2 (g, Q) = {v\\q) 2 + (f±Q) 2 gives the dispersion of the Goldstone mode 8 C , 
and u+(q,Q) = + (fug) 2 + (v±Q) 2 , where A\ = Z^l^^X, the dispersion of 
the longitudinal mode r\ c . Hence we conclude that the BEC phase has two excitation 
branches: a gapless mode, which is expected from Goldstone's theorem, and a gapped 
mode, the longitudinal mode, which describes oscillations of the magnitude of the 
order parameter about its ground state value ip c . It is worth comparing these 
phenomenological results with those of the analysis of gaussian fluctuations (RPA + 
SMA) about the mean- field state described in |Appendix A There we find, in the small 
(q, Q) limit: 

cuUq,Q)^(v ll q) 2 + (vt ) Q) 2 , (40) 

u 2 + (?. Q) * + (v\\q? + ( V ± +) Q) 2 > (41) 

where t»i _) = A ia /k, A 2 + = A 2 (8K - iy 2 /2h 2 , and v^ = A 2 a/h, being A p = 
A s sin(p7r) the breather energy gaps and 



A s = hv s po 



p + 1 fnp\ 4J 
r] tan — 



i+p 



]n 2 (p) \ 2 J v s p 



2(1-P) / J 



PlD 



(42) 



is the soliton gap. One noticeable difference of these results from those of the GL theory 
is that the values of v± predicted by the RPA+SMA for the longitudinal and Goldstone 
modes differ, i.e. ^ v\_ \ This is due to a poor treatment of the Landau- Ginburg 
theory that neglected terms of Cql of order higher than quadratic in the fields r\ c and 9 C . 
These terms represent interactions between the Goldstone and longitudinal field, as well 
as self-interactions of the longitudinal field. On the other hand, the RPA treatment of 
| Appendix A| takes into account some of these interactions and thus yields a correction to 
the dispersion of the Goldstone and longitudinal modes, which leads to different values 
of v±. 

It is interesting to consider the behavior of the system in the limit of strong 
Josephson coupling where J > \1\d- Under these conditions the motion of the atoms is 
not mainly confined to ID and thus they can more easily avoid collisions. Therefore, 
the effect of fluctuations is expected to be less dramatic and assuming the existence of 
a 3D condensate is a good starting point. This leads to Gross-Pitaevskii (GP) theory, 
which, assuming a slowly varying configuration of the order parameter, has the following 
Lagrangian: 

£ GP (x, R, t) = hZ x \W* C (x, R, t)dt* e {x, R, t) - ^-\dj! c {x, R, t) \ 2 



2M 1 



- ^|V R vI/ c (x,R,t)| 2 ] - ^ (|* c (x,R,t)| 2 - |^ C | 2 ) Z .(43) 



2 



nx 
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where Z x = a~ 2 , y± = a 2 J/h, and A = Zig/h, g = 4iih 2 a s J dydz\Wn(y, z)\ A /M 
(Wn(y, z) is the Wannier function at site R). The condensate fraction ip c is the solution 
to the (time-independent) GP equation: 

g\t/j c \ 2 = z c J + fi, (44) 

where \x is the chemical potential. Formally, and ignoring the different values of the 
couplings, Cgp looks almost identical to Cgl except for the fact that ihZi^/*d t ^ c 
replaces Z 2 \d t ^ c \ 2 . However, this change has a dramatic effect on the excitation 
spectrum as we show in what follows. Setting ^ c (x, R, t) = [ip c + i] c (x, R, t)] e^^^t) 
and carrying out the same steps as above, we arrive at: 

C GP (x, R, t) = HZ 1 [- 2iP cVc d t c - ^ (d xVc ) 2 - V -± (V R r/ c ) 2 - 2A^ C 2 

-^(^ c ) 2 -f^ 2 (V^ c ) 2 ]. (45) 

From the above expression, we see that the amplitude r\ c has no dynamics, i.e. there is 
not term involving dtf] c that is not a total derivative ||. The term —2if) c i] c dt6 c , however, 
couples the dynamics of the phase 8 C and the amplitude r] c . This is very different from 
the situation with the GL theory encountered above. We can integrate out r] c , e.g. by 
using the equations of motion for i] c , which yield f] c (x, R, t) = —Zid t 9 c (x, R, t) / (2\tp c )+ 
gradient terms. Thus, the lagrangian becomes: 

£ G p(x, r, t) = S [{d t e c f - vl {d x e c f - v\ (v R # c ) 2 ] , (46) 

being v\ = 7# C | 2 A/ \2MZ X ) and v\ = y ± X\^j c \ 2 /(2Z 1 ). This means that a BEC 
described by the GP theory has a single low-energy mode, whose dispersion, u(q, Q) = 



v 2 q 2 + f^Q 2 , vanishes in the long wave-length limit. This is precisely the Goldstone 
mode, which corresponds to the one that we have already encountered above when 
analyzing the excitation spectrum of the GL theory that describes the extremely 
anisotropic limit where J fim- The presence of the Goldstone mode as the lowest 
energy excitation in both theories is a direct consequence of the spontaneous symmetry 
breaking of global gauge invariance. Thus both kinds of superfluids are adiabatically 
connected as J is increased. On the other hand, the latter theory also exhibits a mode, 
the longitudinal mode, that is gapped for (q, Q) = (0,0). However, the fate of the 
longitudinal mode as J increases is not clear from the above discussion. To gain some 
insight into this issue, let us reconsider the derivation of the GL Lagrangian in the light 
of the previous discussion of GP theory. Apparently, the Lorentz invariance discussed 
above forbids the existence of a term ihZi^/*d t ^ c (x, R, t). Nevertheless, we have to 
remember the emergence of such an invariance in equation ()35)1 is a consequence of the 
linear dispersion of the TLL modes of the uncoupled tubes. The dispersion is linear 
for energies Tiu <^ jiio, but for J ~ ni D non-linear corrections must be taken into 
account. These break the effective Lorentz invariance of (jHB^) and thus a term of the 
|| total derivatives have been dropped from both £qp and Cgl- 
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form i%Zi^f*dt^c is allowed. Phenomenologically, we can consider a modified GL theory 
of the form: 

C' GL (x, R, t) = ihZ^dtVc + [\d t V c (x, R, t)\ 2 - v 2 \d x V c (x, R, t) | 2 

- wilVR^x, R, t)| 2 ] - ™ (|* c (x, R, t)| 2 - ^ff . (47) 

where Z\ is very small for J/h\d <C 1 but grows with J and should grow large for 
J/ Hid > 1- We can perform the same analysis as above to find the excitation spectrum, 
which leads to the following quartic equation for the dispersion of the modes: 



AZ 2 
~Zj 



L0 2 -u±{q,Q<) LU 2 -Lui(q,Q] +^uj 2 = 0. (48) 



J 2 

where uj±(q, Q) are the dispersion of the Goldstone and longitudinal modes for Z\ = 0, 
which have been found above. We see that for small q and Q there is a solution 
such that uj 2 = Aq 2 + BQ 2 . Furthermore, there is also a gapped mode, for which 
us 2 — > A 2 + + (Zi/Z 2 ) 2 as q and Q) tend towards zero. Thus we see that the gap of 
the longitudinal mode depends on the ratio Z 1 /Z 2 . For Z\ C Z 2 , the dispersion of the 
longitudinal mode is essentially the same as that obtained from the GL theory. However, 
as Zx/Z 2 grows with J, the gap increases and therefore the exciting the longitudinal 
mode becomes more costly energetically. At this point is useful to recall that the RPA 
analysis in | Appendix A| indicates the existence of a continuum of soliton-antisoliton 
excitations above a threshold hu = 2A S . The stability of the longitudinal mode is 
ensured if A + < 2A S , which seems to be the case provided one trusts the single- mode 
approximation employed to estimate A + . However, as the gap is pushed up towards 
higher energies by an increasing Zi/Z 2 , the longitudinal mode is likely to acquire a 
large finite lifetime by coupling to this continuum. This will eventually lead to its 
disappearance from the spectrum. 

3.2. Self- consistent harmonic approximation (SCHA) 

An alternative approach to obtain the properties of the BEC phase is the self- 
consistent Harmonic approximation (SCHA) (see e.g. Ref. [T7j). The idea behind this 
approximation is to find an optimal (in the sense to be specified below) quadratic 
Hamiltonian that approximates the actual Hamiltonian of the system. The method 
is more suitably formulated in the path-integral language, where the role of the 
Hamiltonian is played by the euclidean action. The latter is obtained by performing, 
in the Lagrangian of equation ()H5j) . a Wick rotation to imaginary time, t —it, and 
defining the action as S = — Jq dx J^ 13 dT C e s(x,ir) /%. Hence, S[6] = S [6] + Sj[6], 
where 



(49) 



K ^ r L ri 3 
Sole] =^E dx dr 
2n r£ Jo Jo 

Sj{6]= --A B poJ £ / dx / dTCOs(9 R (x,T)-6 w (x,T)). (50) 
n /rra o Jo 



— (d T 0n(x, t)) 2 + v s (d x 9n,(x, t))' 

V s 



Quasi-ID interacting Bose gases 



19 



The partition function reads: 



[dB] e 



In the SCHA, S[8] is approximated by a gaussian (i.e. quadratic) action: 



S v [0] 



G v 1 (qMiU n )\6 Q (q,uj n )\' 



(51) 



(52) 



2M h3L 

where 0q(<7, &> n ) is the Fourier transform of 0r(x, t). This gaussian form is motivated 
by considering the action, equation ([5U]). in the limit of large J. We then expect that 
the replacement cos (0r — 6*r/) ~ 1 — h (0r — #r') 2 would be reasonable. This yields 
a quadratic action and a phonon propagator the form of (|58|). To find the optimal 
G v (g, Q, a> n ), we use Feynman's variational principle 



F = --In Z <F'[G, 

P 



where 



F v - ((S[6] - S y [9]))y, 
.) v stands for the trace over configurations of 9-r(x, t) using e~ Sv ^-' 



and 



-/3F V 



Therefore, by extremizing F'[G V ] with respect to G v , that is, by solving 

5F'[G V ] 
5G v (q, Q,u n 
we find the optimal G v : 



0. 



z c A B (K)p c 



Ja 2 F{Q) e 



2M Lhl3 



(53) 



(54) 



(55) 



(56) 



G Y (q, Q,u n ) G (q,uj n ) 2k 
where F(Q) = St (l — e 1 ^'*), with t the vectors that connect a lattice point to 
its nearest neighbors, and G () 1 (q, uj n ) = u)\ + t> 2 g 2 the free phonon propagator 
(f|l = f s ). This self-consistency equation (|56j) can be solved by rewriting it into two 
equations as follows: defining 

2a 2 



v 



B 



71Z C A 

IK 
we arrive at: 



(Poo) 



Ja\ 
TJ 



G v (q,Q,u n ) 



K 



1TV\\ 



^ n + {v\\qf + vlF{Qt) 



(5* 



The set of equations (|57j) and (|S5|l are solved in |Appendix CD , for the variational 
parameter v\ which is the transverse phonon velocity of the 3D superfluid. 

At T = 0, the equation for Vj_ — v±(T — 0) can be solved analytically. We merely 
state here the result (the details of the calculation can be found in |Appendix Cj ): 

,0 



Ull 



B A B (K)\ (hv\\po\ z c J 



K 



Hid 



Hid 



(59) 



where £>o ~ 1.2971 for a 2D square latttice. Note that the variational propagator 
G y (uj n ,q, Q) has a pole, which for small q and Q, is located at u>^(q, uS) ~ t> 2 g 2 + 
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(i>_i_Q) 2 , for a square 2D transverse lattice. This is the Goldstone mode with v± ~ 
{J/ Hid) 2K ^ AK ~ l \ a power law with an exponent that agrees with the RPA+SMA result 
of | Appendix At where it is shown that 

v ^ ~ Ai ~ A s ~ {J/hid) 2K/(ak - 1] - However, 
the SCHA fails to capture the existence of the longitudinal mode. This is because, as 
described above, the SCHA is tantamount to expanding the cosine in the Josephson 
coupling term. Such an approximation neglects the fact that the phase 9n(x, r) is 2n- 
periodic, and that as a consequence the system supports vortex excitations. On the 
other hand, this is captured by RPA because it does not disregard the periodic aspect 
of the phase since the mean-field state is described by the sine-Gordon model. The 
latter has soliton, anti-soliton, and breather excitations, which are a direct consequence 
of the periodicity of R . 

In the weakly interacting limit (i.e. for K ^> 1), equation (}59|) becomes: 

V l ~ J21, (60) 

V\\ V PlD 

where we have used that Hid / (v\\Po) ~ 7 1//2 = tt / 'K for K ^> 1. Thus, the requirement 
that J / Hid <C 1 automatically implies that v\/v\\ <C 1, that is, we deal with a very 
anisotropic BEC. 

At T > 0, the variational equations must solved numerically (see | Appendix C ) . 



The absence of a solution to equation ()57|) such that for T > T c v ± = allows us to also 
determine the condensation temperature. Thus for weakly interacting Bose gas (K 3> 1 
or 7 small), we find: 

7 ' '" rJ (0.610 + 0.698iT) , (61) 



Hid V f'lD 

This is the same form as that obtained for large K from the RPA analysis of section l3.1.21 
equation (|32|). Thus, the same caveats as those discussed there apply to it. 

At T = T c the velocity v± = has finite a jump, which for K ^> 1 is given by 



Vj_ V V>ID 



(0.376 + 0.274FT x ) . (62) 



This is indeed an artifact of the SCHA, which incorrectly captures the second-order 
character of the transition from the BEC to the "normal" (TLL) phase. 

3. 3. Momentum distribution function 

Besides the excitation spectrum, the SCHA also allows us to obtain the momentum 
distribution of the system in the BEC phase. The latter can be measured in a time-of- 
flight experiment and it is the Fourier transform of the one-body density matrix: 

n(k, K) = jj-j- J dxdx' J dr ± dr' ± e -^±- r '±) e -^-*') ^t( Xj r± )*(x', r' ± )). (63) 

Introducing *$?(x, r±) = H-R^rjJ \E'r(x), where rj_ = (y,z), n(k,K) can be written 
as: 

n(k,K) = £ e -* K '( R - R ') / dxdx'e-^-^ (*l(x)* w (x')) , (64) 

M L R,R' J 
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where VF(K) is the Fourier transform of the Wannier orbital Wo(rj_). Using the 
bosonization formula and within the SCHA 

(¥ R (x)V o (0)) = p A B (K)(e- i[9R{x) - 9om } ~ p ^B(^)e-5 <[eR( "^ e ° {0)l2>SCHA . (65) 

Since i([Ma;)-£o(0)] 2 >scHA = (#o(0)>scha - (0r(x)0 o (O))scha and (M^o(0)> = 
G v (x, R, t — 0), where 



G v (x,R,t) = v^E /^e^ R ^G v (g,Q, 



uj + id) 



(66) 



(67) 



" 8n 2 v 2 ± K (x/v^ 2 + (R/*; ± ) 2 ' 

where G v (q, Q, u n ) is given by equation (|55)l. The result on the last line is the asymptotic 
behavior for large x and R. Thus we see that for large distances G Y (x, R, t) decays very 
rapidly and therefore it is justified to expand e G v(a;,R,t=o) ^ 1 + G v (x, R, t = 0) + . . . 
Keeping only the lowest order term, we obtain the following expression for (k, K) near 
zero momentum: 



n(k,K) ~ \W{K) 



(M L)\i; c \ 2 S kfi S Ki0 + 



v\\Zo 



( Vll k) 2 + (v ± K) 



(68) 



where Z = \ip c \ 2 7i/2K, and the (modulus squared of the) condensate fraction, 



l^cl 2 = PoA B (K) e - G ^=°< R =°) = p Ab(K) 



hv i e 



-Co\V2* 



2fi 1D a 



Po 



(69) 



where Co — 0.2365 for a 2D square lattice. Thus the SCHA approximation also yields 
a condensate fraction that behaves as a power law of J / fiio- The exponent is the same 
as the one obtained in section 13.1.11 from mean-field theory, equation 1)25)1 . 

The momentum distribution can be also obtained at finite temperatures from: 



G v (fc,K,t = 0) 



1 v nvuK- 1 

n^h< + ^\\k) 2 + viF{K) 



ITVnK 1 



coth 



hJ( Vll k) 2 + vlF(K) 



'(v liq y+viF(K) 
Thus for k and K near zero, and T > we find 



T 



n(k,K) ~ \W(K)\' 



(M L)|^(T)|\oVo + 



Tv\\Zo(T)/h 



(70) 
(71) 

(72) 



(v {l k) 2 + (v ± K) 2 _ 

with Zq(T) defined as above with C replaced by ip c {T). 

There are several points that are worth discussion about the results of equations 
and (J72j) . The first is that these form are anisotropic generalizations of the result that 
can be obtained from Bogoliubov theory for a BEC [22] • However, for sufficiently large q 



we expect a crossover to the same power-law behavior of a single ID Bose-gas ~ \k\2K 
(for T = and L — > oo). This is nothing but a reflection of the dimensional crossover 
physics: bosons with kinetic energy e(k) = h 2 k 2 /2M J (but such that e(k) <C H\d for 
the bosonization to remain valid) cannot "feel" the BEC and the intertube coherence. 
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However, they still "feel" the characteristic quantum fluctuations of the TLL state. 
The second observation is related to our discussion of section 13. 1.11 on the experimental 
"coherence fraction" : it is clear by fitting a gaussian to n(k, K) about zero momentum, 
one does not only pick a contribution from \ip c (T)\ 2 but also from the tail of the non- 
condensate fraction. The presence of a trap (and the inhomogeneity of the system caused 
it) should not dramatically modify this conclusion. However, as both contributions are 
positive, we expect that the experimentally measured coherence fraction is indeed an 
upper bound for \ip c (T)\ 2 . 

4. Deconfinement 

In this section we shall obtain the phase diagram of the system for the case where 
the Mott potential is relevant. One necessary condition for this is that 5 = and 
uo 7^ so that g u ^ in equation (JTTJ). As it has been mentioned above, under 
these conditions, and provided repulsion between the bosons is strong enough so that 
the parameter K < 2, the system undergoes a transition to a Mott insulating state 
for J = and arbitrarily small Uq El E]- I n the Mott insulating state 
bosons are localized about the minima of the potential. The question thus is what 
happens as soon as a small Josephson coupling between neighboring tubes is present. 
The latter favors derealization of the bosons leading to the appearance of a BEC, 
and therefore competes with the localization tendencies favored by the Mott potential. 
The competition between these two opposite tendencies leads to a new quantum phase 
transition known as deconfinement. For small values of \uq\ and J, it can be studied 
by means of the renormalization group (RG). This is a much better approach than 
comparing ground state energies of mean-field theories for each phase, as it takes 
into account the fluctuations, which are entirely disregarded within any mean-field 
approximation. However, although the RG approach allows us to estimate the position 
of the boundaries between the Mott insulating and BEC phases, in the present case it 
cannot provide any insight into the nature of the deconfinement transition. For this, we 
shall rely on the mean-field approximation and map the mean field hamiltonian at the 
special value of K = 1/2 to a spin-chain model. 

4-1. Renormalization Group calculation 

Let us return to the low-energy effective field theory defined by the Hamiltonian in 
equation (JTTJ). Since it is a description valid for low-energies and temperatures, the 
couplings K, g u , gj, etc. all depend on the cut-off energy scale, which is typically set by 
the temperature for A = T < p, 1D . Thus, as the temperature is decreased, the couplings 
get renormalized due to virtual transitions to states with energy > A = T, which now 
are excluded from the description. The nature of the ground state is thus determined 
by the dominant coupling as T = A — > 0. 

The change ('flow') for the couplings can be described by a set of differential 
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equations. For small values of g,j and g u the RG flow equations can be obtained 
perturbatively to second order in the couplings. The details of such a calculation are 
given in |Appendix D The result reads: 

dg F gj 



dl K' 

dg 1= ( 9 _±\ n ,9J9F 



(73) 

:j-i:jr Iha\ 

dl V 2KJ 9J + ^K> (?4) 
^ = (2 - K) g u , (75) 

= Ag) - g 2 u K\ (76) 

where I ~ ln/x/T. The coupling gp is generated by the RG and describes an interaction 
between bosons in neighboring tubes. However, since initially there is not such an 
interaction (i.e.g F (0) = 0) and as dg F /d£ = O(gj) its effect on the flow of the other 
couplings is not very important. 

The phase boundary between the BEC and ID Mott insulator phases can be 
obtained by studying the asymptotic behavior of the solutions of equations (fTHJ) to f7H| . 
As the temperature is lowered, i.e. as £ — > +oo, in the regime | < K < 2 % the effective 
couplings g u (£) and gj(£) grow as both the Mott potential and the Josephson coupling 
are relevant perturbations in the RG sense. If g u {£) grows larger, the localization 
tendencies of the Mott insulator phase set in, whereas if gj(£) grows larger, derealization 
tendencies of the BEC phase described in section |3] dominate in the ground state. Both 
phases correspond to two different strong coupling fixed points of the RG. When both 
perturbations are relevant, one has a crossing of the two strong coupling fixed points, 
and since we obtained the above equations assuming that both g u and gj are small, the 
RG is not perturbatively controlled when g u and gj grow large. Note that, given this 
limitation, we cannot entirely exclude that intermediate phases might exist. We will 
come back to that particular point in the next section. 

The relative magnitude of the initial (z.e.bare) values of gj and g u determines 
which coupling grows faster to a value of order one, where the above RG equations 
cease to be valid. The faster growth of one coupling inhibits the other's growth via the 
renormalization of the parameter K: If g u (£) grows faster, K flows towards 0, leading to 
the ID Mott insulator. This reflects the fact that, for large enough g u , the Mott potential 
drives the system towards an incompressible phase, where boson number fluctuations 
are strongly suppressed by the localization of the particles. On the other hand, if gj{£) 
grows faster, K also grows and leads to the 3D BEC phase. A rather crude estimate 
of the boundary can be obtained by ignoring the renormalization of K. In such a case, 
dividing (|74|) by (|75|) . keeping only the leading order terms, and integrating this equation 
up to the scale £*, where both gj(£*) = g u (£*) = s, where s ~ 1 we get that the critical 

% The regime where K < 1 describes a ID system of bosons interacting via long-range forces Thus 
it is not physical for a system of cold atoms interacting via short-range forces. However, the discussion 
in this section will be done for a general value of K. 
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(Y- + ») (y-8) ( Y = 3.5) (y-2) 

K 

Figure 2. Zero-temperature phase diagram of a 2D array of coupled, infinite, ID 
interacting boson systems. The values on the vertical axis are defined up to a factor of 
order unity. K is the Luttinger-liquid parameter (see section l2~2l for explanations). The 
corresponding values of Lieb-Liniger gas parameter 7 = Mg/h 2 po are also displayed. 
A longitudinal commensurate optical potential of amplitude uq (uo ~ 0.03 /iiu for 
K = 1 and ito ~ 0.02 fi for K = 2.5) is assumed. The dashed line corresponds to the 
crude estimate for the phase boundary obtained from equation (|77|l . 



value of g c j = gj(Q) for a given value of <? u (0) scales as: 

2-1/2A' 

0j~0„(O)t=^. (77) 

In the next section we will see that this power-law is also recovered in the mean-field 
treatment of the following section. However, for a more precise estimate of the boundary 
one needs to integrate numerically the RG-flow equations and take into account the 
renormalization of K, which can be important especially near K = 2. This is shown 
in figure where we compare the estimate of (|77|) with the phase boundary obtained 
from a numerical integration of the RG equations. 

The RG equations also allow to extract the gap of the Mott insulator as well 
as the soliton energy gap that characterizes the decay of the one-body correlation 
functions (i.e. the healing length) in the BEC phase. This is achieved by integrating 
the equations until the fastest grown coupling is of order one. One can then use that, at 
this scale £*, the gap A(£) is also of order one, and since gaps renormalize like energies, 
i.e. A(£) = A(£ = 0)e~ e . For example, in the BEC phase, using equation (J72J), keeping 
only the leading term, and integrating this equation up to a scale t ~ log(A/A s ) where 
9j(^*) ~ 1 gives A s ~ [J / I^id) 1 ^ 2 ^ 1 ^ 2 ^-, which has the same power-law behavior as 
the one obtained from the mean-field theory, equation (j42j) . For | < K < 2, a better 
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estimate (numerical) can be obtained by accounting for the renormalization of K using 
the full set of RG equations. However, for K > 2, the renormalization of K is negligible 
and the previous estimate of A s should be accurate. 

Another instance where it can be necessary to take into account the renormalization 
of K is the evaluation of the condensation temperature T c in the presence of the Mott 
potential. In particular, for K < 2, although the initial conditions are such that the 
coupling gj is the dominant one, K may undergo an important renormalization if g u 
is non-zero. The renormalized value of K could be then used in combination with the 
RPA formula, equation (J3Uj) . for T c . 



4-2. Mean field theory of the deconfinement transition 

The RG approach explained above does not provide any insight into the nature of 
the deconfinement transition, e.g. it cannot answer questions such as whether the 
transition is continuous (i.e. second order) or not, or what universality class it belongs 
to. In order to address, at least partially, some of these issues we will treat again the 
interchain coupling in a mean field approximation. We begin by rewriting the mean-field 
Hamiltonian of equation (JT7j) in the presence of the Mott potential: 



rMF 



TT 



K (d x 0(x)f + K- 1 (d x <j>{x)f 



+ 2p uo J dx cos20(x) — 2Jzc\f AbPo\^Pc\ J dx cos 6 (x) + J z c L\ip c \ 2 . (78) 

The new mean-field Hamiltonian is now a double sine-Gordon theory, which is not 
exactly solvable for general values of the parameter K [46J. However, for K = ^ the 
solution takes a particularly simple form as one can directly relate it to the results 
obtained in section 13.11 for the sine-Gordon model. As it has been pointed out above, 
the regime K < 1 is not physical for models like the Lieb-Liniger model or the ID 
Bose-Hubbard model. However, we expect that the solution at this point captures some 
of the essential features of the deconfinement transition for K > 1. Indeed in particular 
in the RG flow nothing special occurs at K = 1 so we can expect the two limits to 
be smoothly connected. The key observation behind the solution at K — | is that, 
as shown in |Appendix E at this point (fTHJ) is the effective low-energy description of a 



Heisenberg anti-ferromagnetic spin chain in the presence of a staggered magnetic field: 
H MF = J £ S m • S m+1 + h • ]T(-ir S m + Jz c L\^ c \\ (79) 

m m 

where h = (h x ,0,h z ), h x = —^/8Trzc(ABPo) 1 ^ 2 \i'c\a'oJ and h z = 27ru (poao), where ao is 
the short-distance cut-off. We can perform a rotation about the OY axis, so that the 
spin variables transform as: 

S x m = cos fS^ + sin <p S z m , (80) 

S z m = - sin cos y?^, (81) 

where tan<£> = h z /h x , and the mean-field Hamiltonian becomes: 

H MF = JoES m • S m+1 + |h| E(-!) m ^ + JzcL\^ c \ 2 . (82) 
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Applying the the spins S m , the bosonization rules the spins given in Appendix E the 
above Hamiltonian becomes: 



H, 



MF 
eff 



2vr 



-{d x e{x)f + 2{d x <P{x)f 



+ 



3) 



f , dx cos 9 (x) + J zcL\ip c \ 2 . 
aov2ir J o 

This is another sG model, similar to the one that we encountered in section 13.11 One 
can again make use of equations (fTTH EQJ) and (f2*Tf to obtain the ground-state energy- 
density, 

S sG (\h\) + z c J\A\ 2 - 



7T.MF 

L 



/L 



The behavior of the order parameter ip c follows from the extremum condition: 
d 



dip c 



S sG (\h\) + z c J\4> c 



0. 



4) 



(85) 



The algebra is essentially very similar to that in section 13.1.11 and we just quote the 
results here. There exists a critical tunneling J c , when the condensate fraction ip c — > 0, 

= 1/2, 

7rp aoKl/3)] 4 



i.e. h x — > 0. We find at K 



ZgJc 
v \\Po, 



2[±77tan(7r/6)F 



'Wo, 



86) 



Thus the condensate fraction grows continuously from zero at the deconfinement 
transition from the Mott Insulator at J < J c to the 3D anisotropic superfluid at J > J c 
according to: 



IPl = POI 



tan 



z c J" 



7) 



Kl/3)] 4 L6 V6/J V^HPo, 
where «(l/3) has been defined in (|22p and at K = 1/2 the dimensionless ratio 
V = PoOo^b(1/2). Note that the scaling J c oc u^ 3 in (jHSj) agrees with that deduced 
from the RG approach of equation (f77)) upon setting K = 1/2. 

It is instructive to compare this deconfinement scenario in 3D with the case where 
only two ID systems are coupled together (via the same Josephson coupling term), i.e. 
a bosonic two-leg ladder, which has been studied in [37] . In both the quasi-lD optical 
lattice and the ladder, there is the competition between the Josephson coupling that 
favors superfluidity, and the localizing tendency of the commensurate periodic potential 
on top of which the interacting bosons hop. This is manifested in a very similar structure 
of the bosonized Hamiltonians. For the ladder case, it is convenient use a symmetric 
and anti-symmetric combination of the fields of the two chains, which leads to slightly 
different forms for the RG equations. Nevertheless, the mentioned competition occurs 
for the same range of K: 1/4 < K < 2, where, in the ladder case, K characterizes to the 
correlations of antisymmetric fields. Qualitatively, the critical intertube coupling for the 
deconfinement transition has the same power-law dependence on g u as in equation (|77jl . 
However, as the ladder is still effectively ID, the nature of the deconfinement transition 
is expected to be very different from the one studied here, where a fully 3D (anisotropic) 
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BEC results. In Ref . jU] , the deconfinement transition in the ladder was found to be in 
the Berezinskii-Kosterlitz-Thouless universality class, at least when (for fixed interaction 
strength U) the intertube hopping is larger than the intratube hopping (note that this 
limit has not been studied in this work as it would correspond to coupled 2D systems, 
with very different physics compared to the coupled ID systems.) 

Note that in the entire phase diagram (in the thermodynamic limit, see figure |2J), 
we have found evidence for the BEC and the ID Mott insulating phases only. We 
have not found any evidence (or hints of evidence) for more exotic states such like 
the sliding Luttinger liquid phase jlHl EH ED] or some kind of supersolid phase. The 
sliding Luttinger liquid [JE1 EEE EH! * s a perfect insulator in the transverse direction (like 
our 2D Mott insulator for finite size systems); along the tube direction, the system 
is still a Tomonaga-Luttinger liquid, but the parameters K and v s are now functions 
that depend on the transverse momentum. This state may be stabilized by a strong 
interaction of the long- wave length part of density fluctuations between tubes [!§]• We 
note that, although the RG procedure described above generates such a density-density 
interaction between the tubes, in the bare Hamiltonian there are no such interactions 
to start with (atoms in different tubes have negligible interactions) and therefore, the 
coupling cjp is likely to remain small as the RG flow proceeds to lower energies, at least 
within weak coupling perturbative RG we use here. Thus it seems unlikely that the 
system can enter the regime of a sliding Luttinger liquid EH] phase. An even more 
exotic possibility would be a phase where the tubes are ID Mott insulators (in the 
presence of a commensurate periodic potential and strong enough interaction in-tube), 
while at the same time, there would be coherence between the tubes in the transverse 
direction. In this scenario, presumably both the couplings g u (Mott potential) and J 
(Josephson coupling) would be relevant. We cannot rule out such a possibility in the 
strong coupling limit. However, we have no evidence for it in our weak coupling analysis. 
Nevertheless, in a quasi- ID system and on physical grounds, it seems hard to imagine 
a phase where the boson density is kept commensurate within each tube, as required 
by the existence of the Mott insulator, while at the same time bosons are able to freely 
hop to establish phase coherence between tubes (but not within tubes!). 

As far as the competition between the Josephson coupling and the Mott potential 
is concerned, mean field theory is a useful tool (in higher dimensions as well as when 
the tubes are coupled) for providing insights into which possible phases and what 
their properties are. However, it is usually inappropriate for predicting the correct 
critical behavior. For cold atomic systems, because the systems are finite and of the 
existence of the harmonic trap that makes the samples inhomogeneous, it is not very 
realistic to study critical behavior as the latter is usually strongly perturbed jSlj. Mean 
field theory seems thus perfectly adapted to the present study. Let us however more 
generally comment on the universality at the transition. Because the field cos(20r) 
is a vortex (instanton) creation operator J7j for the field the Hamiltonian for 
each tube R can be faithfully represented by a classical 2D XY Hamiltonian involving 
J X>T cos[6'r,(x / , t') — 9-r(x,t)]. It thus becomes apparent that the coupling between the 
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tubes leads to an anisotropic version of the 3 + 1 dimensional XY model. The transition 
will thus be in the universality class of the classical D = d + z = 4 XY model in 
agreement with general considerations on superffuid to insulator transitions [39J. In the 
present case, the transition is tuned by varying the anisotropy of the model 5 ~ J — J c . 
For such a transition, the superfluid density of the classical model is expected to behave 
as S u ^ d+Z ~ 2 ', where for the commensurate Mott insulator to superfluid transition z = 1 
and, at d — 3, which is the upper critical dimension, the critical exponent v — 1/2 
is mean-field-like. One would thus find that the superfluid density behaves as 5, in 
agreement with our result for K = 1/2 of equation (f57|). upon identifying the superfluid 
density of the classical XY model with the square BEC fraction of the quantum model, 

Note that in the absence of the self-consistency condition ()18)) . the mean field 
Hamiltonian is equivalent to an XY model in presence of a symmetry breaking 
field jSH E2] • It has been recently argued that such a model can have a sequence of two 
phase transitions |53*| I54j . Since in the present case the coefficient of the Mott potential 
(~ cos(20r)) term is small, the system would be in the regime where it undergoes 
a single phase transition even for the model with the symmetry breaking field. But, 
more importantly, the model that we have studied has to be supplemented by the self- 
consistency condition, equation (fTSj) . which stems from our mean-field treatment of the 
original Josephson coupling cos(6*r — 9r>). This condition will modify the properties of 
the phase where this coupling is relevant, and for the reasons explained above, should 
bring back the model in the universality class of the 3 + 1-dimensional XY model. 

5. Finite size effects 

Actual experimental systems consists of finite 2D arrays of finite ID systems (tubes). 
The finite size of the sample has some well known consequences in condensed matter 
physics such as inducing the quantization of the excitation modes and turning phase 
transitions into more or less sharp crossovers. The latter can be understood using the 
renormalizat ion-group analysis of section 14.11 The finite size of the tubes introduces 
a length scale into the problem, namely the size of the tube L, which effectively cut- 
offs the RG flow at a length scale of the order of L (or, for that matter, the smallest 
length scale characterizing the size of the sample). Thus, flows cannot always proceed to 
strong coupling where the system acquires the properties of one of the thermodynamic 
phases discussed above. Furthermore, as we show below, the finite size L leads to the 
appearance of a new phase. 

5.1. Finite size tubes and 2D Mott transition 

Typically in the experiments [01 E| a number Nq of bosons ranging from a few a tens to a 
few hundreds are confined in a tube of ~ 10/im long. This leads to energy quantization of 
the sound modes in longitudinal direction, but also to a finite energy cost for adding (or 
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Figure 3. Zero-temperature phase diagram of a 2D array of coupled, finite, ID 
interacting boson systems. The values on the vertical axis are defined up to a factor of 
order unity. K is the Luttinger-liquid parameter (see section l2~2l for explanations). The 
corresponding values of Lieb-Liniger gas parameter 7 = Mg/h 2 po are also displayed. 
A longitudinal commensurate optical potential of amplitude uq (uq ss 0.03 hid for 
K = 1 and ito ~ 0.02 /1 for K = 2.5) is assumed. The vertical dashed line describes 
the boundary between the (finite-sized) Tomonaga-Luttinger phase. 



removing) one atom to the tube. This energy scale is Eq = fmv s /KL [21EZ1- The tube 
can be thus regarded as a small atomic "quantum dot" with a characteristic "charging 
energy" E c . If E c is sufficiently large, it can suppress hopping in an analogous way as 
the phenomenon of Coulomb blockade suppressing tunneling through electron quantum 
dots and other mesoscopic systems. In order to know when this will happen Ec must 
be balanced against the hopping energy Ej. The latter is not just JN as one would 
naively expect for non-interacting bosons. Interactions and the fluctuations that they 
induce in the longitudinal direction alter the dependence of Ej on N as we show in 
|Appendix F 



Let us consider a system of finite tubes in the BEC phase (i.e. the Mott potential 
is irrelevant or absent). In the limit of small J (J/fi\D *C 1) of interest here, the 
longitudinal fluctuations with wavenumber q 7^ can be integrated out as discussed in 
Appendix F The resulting effective Hamiltonian is a quantum-phase model identical 



to the one used to describe a 2D array of Josephson junctions: 
H QP = - Ej cos (^or - 0orO 

(R,R'> 

+ t^E(^-^o) 2 -^E^ 

R R 
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where Nr is the particle-number operator of tube at site R, and #or the canonically 
conjugate phase operator: [0 O r, N-r/] = z<5r,r'. The renormalized hopping Ej = 
Ej(N ) ~ JN 2K . The model (jHHj) has been extensively studied in the literature 
(see e.g. Ref. [22] arid references therein). It exhibits a quantum phase transition 
transition between a two-dimensional superfluid (2D SF) phase (a BEC at T = 0) 
and a two-dimensional Mott insulator (2D MI) at a critical value of the ratio Ej/Ec- 
For commensurate filling (integer Nq) using Monte Carlo (MC) van Otterlo et al. [HE] 
found (Ej/Ec)c — 0.15. Using the forms for Eq and Ej derived above, and assuming 
that the tubes are described by the Lieb-Lininger model, the MC result reduces to 
Jc/l^iD — 0.3iV 3 ^ 2 in the Tonks limit (K = 1) and to J c / [1\b — 0.15Nq 2 for weakly 
interacting bosons (K ^> 1) [20] ■ A similar result has been more recently obtained in 
Ref. [Sn] using the random phase approximation (RPA). Our result, however, includes 
critical fluctuations beyond the RPA as it relies on the numerically exact results from 
MC simulations. In figure E] we show the complete zero-temperature phase diagram 
including the possibility of a phase where the tubes are effectively decoupled in a 2D 
MI phase (labeled Decoupled tubes in the diagram). Note that in the 2D MI, only 
the phase coherence between different tubes of the 2D lattice is lost. However, apart 
from the finite size-gap ~ hirv s /L, there is not gap for excitations in the longitudinal 
direction. This makes this phase different from the ID MI described above, where there 
is a gap (depending very weakly on the size L) for longitudinal excitations. If the Mott 
potential is applied to the system in this phase, it will cross over to a ID Mott insulator 
for a finite, but small value, of the dimensionless strength g u . This crossover should 
take place for K < 2 [21 EEl EE3 IH] and is thus represented by a vertical dashed line in 
figure El 

In the superfluid phase, the lattice will also exhibit a crossover as a function of 
temperature: If the condensation temperature T c found in section 13.1.21 is larger than 
the finite-size gap hv s /L, the gas should behave as a 3D BEC and cross over to a 
2D SF for T ~ hv s /L. However, if T c <C hv s /L the system should behave as a 2D 
quasi-condensate with no long range phase coherence (and therefore ip c = 0) at finite 
temperature. Since Tiv s /L ~ Tiojq, u q being the longitudinal trapping frequency, and 
experimentally p[] ujq ~ 0.1 kHz whereas T c /h ~ 10 kHz, we expect to be in a regime 
where the BEC is always observed. The presence of harmonic confinement means that 
at least some of the tubes will have a filling deviating from commensurate filling. The 
critical value of J / fiir, is reduced relative to the half-filled shown in figure HI 

Finally, let us estimate the typical value of the critical ratio J c / Hid- Using the 
above results and assuming that [1\d ~ 0.1 hu± [7\ and N = 100, we conclude 
that J c /huj± ~ 10~ 4 for K — 1 whereas J c /Ttio c ~ 10~ 5 for K ^> 1. The deepest 
lattices used in the experiments of Refs. [3 HO] had V± = 30 Er, which corresponds 
to J/hu± = 5 x 10^ 5 . This is of the same order of magnitude as J c /hujj_ for K ^> 1, 
although a larger number iVo will further reduce this estimate, and therefore for this 
lattice the system may be well in the BEC phase (see figure E])- However, in the 
experiments [3 HO] no coherence was apparently observed between the tubes. The 
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Figure 4. Phase diagram of a 2D array of finite ID Bose gases as a function for the 
local chemical potential fi and the Josephson coupling J for different values of the 
Luttinger parameter: K = 1 (7 — > +00, continuous curve), K = 2 (7 = 3.5, dotted 
curve), = 4 ( 7 = 0.71, dash-dotted curve). 



explanation for this contradiction is that the value of tunneling time h/J must be 
compared with the duration of the experiment t exp and it turns out that for a lattice 
depth of Voj_ = 30 E R this time is indeed longer than t exp . Therefore, over the duration 
of the experiment almost no tunneling process occurs and the build-up of intertube 
coherence can thus hardly take place. 



6. Comparison to experiments 



In this section we discuss and summarize our results in connection with the recent 
experiments on anisotropic optical lattices [2U 13 CHI 1221 UH US] • The complete phase 
diagram, taking into account the finiteness of the ID gas tubes, is shown in figure 01 At 
zero temperature, we predict three phases: an anisotropic BEC, a ID Mott insulator, 
and a phase where the gas tubes are decoupled (i.e. phase incoherent), namely a 2D 
Mott insulator. 

However, as mentioned in section [21 the experimental systems are indeed 
harmonically confined. Unfortunately, so far it has proven very difficult to deal 
analytically and simultaneously with strong correlations and the harmonic trap. On 
the other hand, numerical methods can tackle this situation [SOI EH E21 ESI EH ESj- 
In the harmonic trap, the cloud becomes inhomogeneous, and thus as demonstrated 
in a number of numerical simulations |3S1 EH EH1 EH ESI, f° r a wide range of system's 
parameters, several phases can coexist in the same trap. These results have confirmed 
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the experimental observations ^2] as well as previous theoretical expectations based on 
mean-field calculations [21]. The numerical simulations have been mostly restricted to 
ID systems |HH EHj an d isotropic lattices [M], we expect most of their conclusions to 
also apply to the anisotropic optical lattices that interest us here. Thus, we also expect 
coexistence of the phases described above when the system is confined in a harmonic 
trap. Qualitatively, this indeed can be understood using the local density approximation, 
assuming that the trap potential gives rise to a local chemical potential which varies 
from point to point. Thus, from strong interactions we expect the formation of domains 
of the ID Mott insulating phase near the center of the trap, whereas the superfluid 
phase should be located in the outer part of the cloud. If the transverse lattice is 
very deep (i.e. if Vq± is large) there also can be regions where the number of boson 
per tube is small and the the phase coherence is lost because the charging energy Ec 
becomes locally larger than Ej. Nevertheless, confirmation of this qualitative picture 
will require further experimental and numerical investigation, which we hope this work 
helps to motivate. The coexistence of different phases may be revealed by analyzing 
the visibility of the interference pattern observed in an expansion experiment |23 EE]- 
Let us finally mention that the existence of the trap can also lead to new and poorly 
understood phenomena such as suppression of quantum criticality [51] . 

Taking into account the above, we only aim at a qualitative description of some 
experimental observations. Thus, we have shown in 13.1.11 that the strong quantum 
fluctuations characteristic of ID systems strongly deplete the condensate fraction in 
the BEC phase that forms for arbitrarily small intertube Josephson coupling. The 
condensate fraction can be as low as ~ 10 % for typical experimental parameters (Voj_ = 
20 Er). We have also argued that the experimentally measured coherence fractions 
may also include some of the non-condensate fraction, and therefore, provided heating 
effects of the sample in the measurement can be disregarded, the coherence fraction 
should be an upper bound to the BEC fraction. Quantum and thermal fluctuations are 
also responsible for the power- law behavior of the condensation temperature with J / il\b- 
Estimates of the latter, when compared to the experimental estimates of the temperature 
of the lattice [7] and the finite-size excitation gap means that for Vo± = 20 En the 
system should be in the BEC phase described in section El provided the Mott potential 
is sufficiently weak. In this phase, the system exhibits two kinds of low energy modes in 
the long wave-length limit: a gapless Goldstone mode, and a gapped longitudinal mode. 
Whereas the former corresponds to fluctuations of the phase of the order parameter, 
the latter is associated with fluctuations of its amplitude, namely the BEC density. 

If the strength of the Mott potential is increased, the system (or at least part of it 
near the center of the trap) should undergo a deconfinement quantum phase transition 
(or crossover, given the finiteness of the sample) to a ID Mott insulating phase. This 
phase is characterized by the existence of a gap in the excitation spectrum which is 
largely independent of the size of the tube. The hopping between the tubes is suppressed 
for temperatures (or energies) below this gap, and it is incoherent above it. For a deeper 
optical lattices, we predict that hopping between tubes can also be suppressed by the 
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"charging energy" of the (finite-sized) tubes. We also notice that in some experiments 
the duration time of the experiment is also a limiting factor for the establishment of 
phase coherence even if hopping is not effectively suppressed by the charging energy. 

Stoferle et al. [TU] have studied an optical lattice system for several values of the 
depth of the transverse optical lattice potential (hence J). By analyzing the width 
of the momentum peak near k = in the interference pattern after free expansion, 
they were able to deduce the energy absorption in response to a time-dependent 
modulation of the Mott potential. From this analysis, they could tell when the system 
behaves as a superfluid or as a Mott insulator. These experiments were analyzed 
recently |4"2*| EOl EU EE] and we will not repeat the conclusions of the analysis here. 
In the same series of experiments, they also studied the dependence of the width of the 
momentum peak near k = in the ground state on the strength of the Mott potential, 
for several values of Vo±. For a BEC or a ID superfluid, the inverse of the width is set 
by a geometrical factor and it is proportional to the size of the size of the cloud (or, in 
the ID case, the thermal length hv s /T, whichever the shortest). For a Mott insulator, 
however, the inverse of the width is set by the correlation length £ c ~ hv s /A c , where 
A c is the Mott gap. Thus, an increase of the width as a function of the Mott potential 
reveal the opening of a gap in the spectrum of the system [33] . The experimental curves 
of the width as a function of the Mott potential present upturns around the points where 
the system should enter the Mott insulating phase in the thermodynamic limit. Our 
expectation, based on the phase diagram of figure El is that the deconfinement phase 
transition should take place for values of the Mott potential slightly larger than in the 
pure ID case (represented by the vertical line at K = 2 in figure EJ). It is necessary to 
recall that, in the Bose-Hubbard model (cf. equation (JBJ)) that describes the systems 
of Ref. jni, the Mott potential controls not only the dimensionless parameter g u that 
enters in the RG equations of section 14.11 but also the ratio of the interaction U to the 
kinetic energy J x , and therefore the Luttinger parameter K as well. The latter decreases 
towards the Tonks limit (K = 1) with increasing U/J x . The theoretical expectation is 
in good qualitative agreement with the experimental observation. 

The finite extent of trapped cloud (either longitudinally or transversally) limits 
the minimum momenta of the modes. Thus the energy of the lowest modes in the 3D 
SF phase can be directly estimated from our results by using the minimum available 
momentum in (|4U|41|) . For instance, for a finite 2D lattice containing M y x M z tubes 
{%. e. an atom cloud of size L x M y a x M z a) the lowest available momentum is ~ n /(Mia). 
Putting this value into (j4UI41j) shows that the frequency of the lowest transverse modes 
decreases with decreasing J. This analysis neglects the possibility of phase coexistence 
described above, which calls for further investigation of these issues. 

We have been not able to calculate the finite temperature excitation spectrum 
across the dimensional crossover, but it is interesting to speculate what happens 
to the Tomonaga-Luttinger liquid (TLL) spectrum, especially when the longitudinal 
momentum q is near 27rpo- On general grounds, at temperatures large compared to 
the dimensional crossover scale ~ T c , the system should exhibit essentially the ID 




q 

Figure 5. Schematic excitation spectrum of the BEC phase. Shaded lines and curves 
represent different types of excitations. The frequency of the excitation is ui and q is 
momentum in the longitudinal direction {i.e. parallel to the tube axis). The continuum 
of particle and hole excitations stems from the spectrum a ID system, as found by 
Lieb in Ref. [SSI- At lower energies, however, the system develops the phase coherence 
characteristic of a BEC and therefore must obey Landau's criterion of superfluidity. 
In particular, this means that the system must exhibit a gap at q = 2irpQ, where po 
is the linear density of the ID system. In the long wave-length limit, we find the 
system exhibits two kinds of excitations: a gapless (Goldstone) mode, and a gapped 
longitudinal mode. There is also a continuum of excitations above the longitudinal 
mode (see discussion at the end of section 

spectrum of a TLL, with its characteristic continuum of (particle-hole) excitations 
extending roughly between q « and q ~ 2np [63J. As the temperature or the 
excitation frequency get below the crossover scale, the low-energy excitations around 
2np must disappear: being the BEC a 3D superfluid, by Landau's criterion, there 
can be no low energy excitation that contributes to dissipation of superfluid motion. 
The gap about q = 27rp in the spectrum of the anisotropic BEC phase would be thus 
reminiscent of the roton gap in liquid Helium, which may be regarded as the remnant of 
the crystallization tendency at the momentum of where the roton minimum occurs [64|. 
In our case, because of the underlying ID physics, the roton gap should have a continuum 
of excitations above it, which is the remnant of the ID continuum of excitation in the 
TLL (see figure EJ). This may account for the "relatively broad" energy absorption 
spectrum observed in the "ID to 3D crossover" regime in the experiment of Stoferle et 
al. PHI12] 

7. Conclusions and discussion 

In summary, through a variety of methods, we have arrived at the following picture 
For the case where the Mott potential is either absent, incommensurate, or generally, 
irrelevant in the RG , the existence of an arbitrarily weak intertube tunneling 
immediately turns the ID arrays into a 3D BEC at zero temperature, albeit with 
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anisotropic properties. At finite temperature, thermal fluctuations cause a phase 
transition at T = T c where the 3D coherence (and the BEC) is lost. Since T c <C {i\D 
in the quasi-lD optical lattices studied here, for T c < T ^ hid the system behaves 
as a 2D array of phase-incoherent Tomonaga-Luttinger liquids. We have obtained 
explicit expressions for T c , the condensate fraction at T = 0, the excitation spectrum, 
and momentum distribution function, which could be compared to measurements in 
the experimental systems. In particular, the power-law behaviors of J / fi found for 
T c and condensate fraction are consequences of the strong fluctuations that dominate 
the properties of the uncoupled ID systems. In the presence of a relevant (in the RG 
sense) Mott potential, a finite intertube tunneling (J c ) is required to overcome the 
localization tendencies of the interacting bosons. This process takes place in the form 
of a deconfinement quantum phase transition. Our main result in this case is the phase 
diagram shown in figure 01 as a function of J/ \i\r> and the parameter K characterizing 
the interactions in the tube. For the deconfinement transition, we also find a power-law 
dependence of J c on the strength of the Mott potential Uq, and a power-law dependence 
of the condensate fraction near the transition (indeed, these two results are strictly 
valid for K = 1/2, which is outside the range of values for bosons with short-range 
interactions, but we speculate that these results are adiabatically connected to the 
physical (for bosons) regime of K > 1 as we have found no evidence of any quantum 
phase transitions for 1/2 > K > 1 in our RG calculation.) Finally, for comparisons 
to experiments, we have looked into trap effects by studying tubes of finite length and 
described the transition where the system goes from a system of uncoupled ID Luttinger 
liquids to the 3D superfluid as (renormalized) J overcomes the "charging energy" of the 
finite length tube. As discussed in details in the previous section we can compare our 
predictions with the experiments of Stoferle et al. [10J. Some of our predictions, in 
particular the large reduction of the coherent fraction is in qualitative agreement with 
the experiment. 

Clearly, further experimental and theoretical studies are needed to provide a more 
thorough and quantitative picture of the dimensional crossover and the deconfinement 
transition. From the theoretical point of view, taking into account the harmonic trap in 
a more quantitative way would be a step towards a more quantitative description of the 
experiments. From the experimental point of view, investigation of larger systems with 
more atoms per tube and more tubes, and in longer experiments that allow to probe 
smaller values of the Josephson coupling J, would permit a more complete exploration 
of the phase diagram of the quasi-lD lattices. Also, most of the current experiments in 
these lattices have been carried out using rather deep Mott potentials in order to reduce 
the ratio of kinetic to interaction energy (this is the regime where the Bose-Hubbard 
model of equation (JBJ) is applicable). However, It would be very interesting to conduct 
experiments where the interaction energy between the atoms is tuned independently of 
the longitudinally applied (Mott) potential. Increasing the interaction strength in the 
absence of Mott potential was recently achieved by Kinoshita et al. by confining 
the atoms in tighter ID traps. Under these conditions, and for interaction strength 
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7 ~ 3.5 (21 E] in ID, the application of the Mott potential of strength much smaller 
than the chemical potential \L\b should suffice to cause localization of the bosons in 
the ID Mott insulating phase. To the best of our knowledge, this "weak lattice" Mott 
insultator has not yet really been realized. In this phase quantum effects are stronger, 
and it would be then very interesting to drive the system across the deconfmement 
transition by changing J. 
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Appendix A. Mean-field theory (MFT) and gaussian fluctuations (RPA) 

In this section we shall employ the functional integral to derive the mean-field 
Hamiltonian theory of section 13. 1\ and to go further by taking into account gaussian 
fluctuations about the mean-field state. The method is known as the random-phase 
approximation (RPA). Our treatment follows Ref. |65j . 

Let us consider the model defined by equation (JTTJ). Using the coherent-state path- 
integral, the partition function can be written as follows: 



T~t(i[>i&, V>r) being the Hamiltonian density corresponding to equation (fTTj). Since the 
Mott potential plays no role in the BEC phase (i.e. it is irrelevant in the RG sense), we 
drop it from the Hamiltonian 7i(V , R, ; V'r)- However, the action S[tp*,ip] still contains 
the non-trivial Josephson coupling between neighboring ID systems. 




(A.l) 



where the (euclidean) action, 





(R,R'> 



(A.3) 



R,R' 




(A.4) 



In the above expression we have introduced the following notation: 




(A.5) 



t 
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where t runs over set of vectors that connect a lattice point to its nearest neighbors. 
We perform a Hubbard-Stratonovich (HS) decoupling of Sj [ip* , ip]: 

To make sense of Jr^r/ the Fourier transform of Jr 5 r/ is needed, 

J(Q) = JJ2 e iQ t = 2 J [cos Q y a + cos Q z a] . (A.7) 
t 

The last result applies to a square lattice (in the OYZ plane) with lattice parameter a. 
Hence, 

I e iQ.(R-R') 

J R TV = TF H 777^ • ( A - 8 ) 

Thus, after performing the HS decoupling, the action becomes: 

rh/3 dr r L 

S[?,£,1,;1,] = Sb[PM+Sj[C,Z]+ x / ^E[^^r)+c.c.],(A.9) 

Jo h Jo R 



where 

rh/3 dr r L 



R L 



^2 

^R (^r - A*) ^R + I^Rl 2 + ||^r| 



(A.10) 



5j[r^]=/ T ^E^,r)J R i R 4'(v). (A.ll) 
Jo n Jo 

To perform the mean field approximation we formally integrate out the boson degrees 
of freedom described by ipn,(x,r) and ?/> R (x, t), and obtain the effective action for the 
HS fields: 

SQL[?,t] =Sj[C^] + S^[C^}, (A.12) 

S^Cit] = -ln(e-/o W x/o i& fe^ T )^^ T ) + ^^ T )^M)^^. (A.13) 

We have introduced the notation S GL for the effective action to indicate that it is indeed 
the GL functional. 

To obtain the mean-field equations we perform a saddle-point expansion. The 
saddle-point is determined by the following equations: 

" ft 1? Jr - r ' eR ' c( "' r) + ^R(^y " °' (A - 14) 

*(R,s,r) - ft^^(^r)J R ,_ R + ^ (XjT) -0. (A.15) 

We look for uniform solutions to the above equations, i.e. £r(x, t) = £ c and £ R (x, t) = 
£*. Furthermore, because of global gauge invariance, the particular choice for the phase 
of £o should not matter, that is, if (£*,£ c ) is a saddle point then so is (e - * 6 ^*, e l9c ^ c ). 
Thus [C?^c] = S^(\C, C \ 2 ). This implies that the above equations reduce to: 

1 ^GL(ie c | 2 ) = lE^R+ir74l-7i^(^ 2 )= ' (A.16) 



M Lh{3 1 ' h V 
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where the prime stands for the derivative with respect to = |£| 2 . Furthermore, using 
equation (jA.13|) . equation (jA.14|) can also be written as: 

£ ^R-R'£ C = -(M*,r)) = -A, (A.17) 

R' 

where it is understood that the expectation value of ipu(x, r) is taken by tracing out 
the boson degrees of freedom in the presence of the constant HS (£*, £ c ), i.e. . using 
•S'tC) £cj V 1 *) (cf- equation (|A.9|1 ). Therefore, ()A.15|) must be solved self-consistently 
for (Cc)Cc) : This is the self-consistency condition mentioned in section l3~Tl Upon 
multiplying equation (|A.17|) by Jr» r from the left and summing over R, we arrive 
at: 

£ C (R, x, t) = £ c = -zctyc, (A. 18) 

where zc is the number of nearest neighbors (zc = 4 for a 2D square lattice). The above 
system of equations defines the mean field theory, which we have solved, in Hamiltonian 
form, in section I3~T1 

We next take into account fluctuations by considering the gaussian corrections to 
the saddle point (£*, £ c )- Let 

5£ R (x,T)=£ R 0r,-r)-£ c , (A. 19) 

8&&t)=&.(x,t)-C (A.20) 

To proceed further we need to integrate out the boson degrees of freedom and obtain 
ity [£*,£] explicitly. Since we are interested in the properties at long wave-lengths and 
low frequencies and temperatures, we first apply bosonization and obtain: 

S[C,t,0\ = J2 S n[0}+SAC^]+^/ABPo / -r / &(*, r)e**<^ + c.c. .(A.21) 

.70 lb J 

In the above expression, 

K r % P , /- L 
2?r Jo 



— (MrO, ^r)) 2 + w s {d x 9^{x, r)) 2 



(A.22) 



In this (bosonized) form, integrating out the boson degrees of freedom amounts to 
tracing out the phase field #r(x, r). This will carried out in what follows using the 
cumulant expansion: 

ln ( e ~ Ar '^L = ^ A [C^,0])sG + \({A [C^B] - (A [C\^0}) sG ) 2 ) sG + ■■■ (A.23) 

where (. . .) s q denotes the trace over (smooth) configurations of 6 l R (x, r) using the weight 
e -s[£*,£ c ,0] e q Ua tion ( IA.21jl . note that this action defines a collection of independent 
sine-Gordon models, which explains the notation). The RPA corresponds to keeping 
the terms of the above expansion up to quadratic order in 5£ and 51;*. Thus, using that 

A[C,^9} = j o -j o dx^[C R (x,r)e w ^ + e- w ^Ux,r)}, (A.24) 
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Po^i 



h/3 



dr 



' dx £ [5&(x, t) (e ie ^) + (e~ w ^) S^x, r) 



R 



fTi/3 P Li 

■ drdr' dxdx' \{{e- ie ^ x ^e- ie ^ x '^))S^(x,r)5Cn'(x', 

J ° J ° RR' L 

+ 2{(e ie ^e- ie ^ x '^)) 5&(x, r) Sfatf, r')] +■■■ 
We have introduced the following notation: 

((F[e(x,T)G[e(x',T f )])) = ((F[6(x,r)] - (F[0(O, 0)]> sG ) 

x (G[9(x',r')]-(G[e(0,0)])s G )) S G, 



2h 2 



(A. 



with 



and 



J[d6] F[0{x,T)]e- s ^ e && 
J[d6] e-^Gld&M 



(A.26) 
(A.27) 



' ' ' dr 



2tt Jo 



+ V A B po 



dxJ2 

R 

o h Jo 



(d T 9 R (x,r)) 2 + v s (d x e K (x,r)f 



dxJ2 [C e WRix ' T) + c.c. 



R 



(A.28) 



Since in this action the different lattice sites are uncoupled, the correlation functions 
are diagonal in R and R'. Thus we can write the quadratic action that describes the 
(gaussian) fluctuations around the saddle point as follows: 

1 r h P r r x 

S ™ [E ' E] = 2hl dT dT ' J dX dX ' ^, \- Er{X > t)3 *-^ X - X '' T ~ t '^ x '> r ') 

+ (s R (z, r) + Ql(x, r)) G R R ,(x -x',r- r') (Z R '(x', r') + Q R >(x', r')) ] (A.29) 

Note that the terms in Sq^ a linear in S and vanish because we expand around the 
saddle point. We have also dropped the constant term that is proportional to the mean 
field free energy. In the above expression the following spinors have been introduced: 



Hr(x,t) 
2 r (x,t) 



,Qr(x,t) 



$£r(x,t) 



Qr(x, t) 
q^(x, r) 



(A.30) 

(A.31) 
(A.32) 



Qb.{x,t) qn{x,T) 

where we have introduced the sources q^{x, r) and gR, x, r), which in the original action 
are coupled to ipR,(x, r) — ip c and ^ R (x, r) — ip*, respectively. The matrices: 

Jr-r'(1 " 2) = J£. w 5(x - x')S(r - r') diag [1, 1] , (A.33) 



G r _r/(1,2) 



(A.34) 



Quasi-ID interacting Bose gases 



40 



with 1 = ( R, x, t) and 2 



^r(i,2) 



R', x 1 , t 1 ), and (p, q 



AbA) //^ R (x J r) (0 i g e R (a:',T') 



ft 



))sG ^R-R',0 



(A.35) 

-0c - 



To perform the averages using a gauge transformation is applied #r(x, r) 
^r(^) r ) so that £ c = —zqJ^c can be treated as a real number. 

Let us now pause to consider the consequences of global gauge invariance on the 
correlation functions of (jA.35|) . As we have found above, the GL functional is a function 
of |£ c | 2 only. Thus, if we expand £V(|£| 2 ) around the saddle point solution: 



OS, 



OS, 



1 d 2 S„ 



2! di 



+ 



1 d 2 S 4 



-Mc\ z )\m z + 



where 



d 2 S 4 
d£ 2 
d 2 S i> 



sum 2 + i&raie 



(A.36) 

(A.37) 

(A.38) 
Thus, 



Using equation ^M, ^(1^1) = -(M„L/3) Er/r = -(M L/3) [J(Q = 0)] 
by setting <5£*, <5£ = const, and taking (3 — > oo in equation ()A.25|) . and comparing with 
equation (jA.36|) . the following results can be obtained: 



Uq = 0,w = 0) = = 0,w = 0) 

9i(q = 0,w = 0) = = 0,w = 0) 

where we have used that lim^ 



<9 2 ,S, 



1 d 2 S,. 



(Kcl 2 ) 
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sG 



n 



d 2 s 



sG 



(A.39) 
(A.40) 



M L(3 d£d? 

S4\£ c \ 2 )/(M Lh{3) = £ sG (|£d 2 ) is the energy density 
of the sine-Gordon model, equation (JI5j). We shall use the above identities to calculate 
of the energy dispersion of the modes in the BEC phase. This phase exists below 
the critical temperature (T c to be determined below), and is characterized by having 
i>c oc £,c 7^ so that global gauge symmetry is spontaneously broken. To compute the 
excitation spectrum in this phase, we diagonalize the matrix J -1 + G (henceforth we 
drop the index R — R'). J -1 is already diagonal, whereas the matrix G is rendered 
diagonal by: 



U 



1 

72 



-iB c 



— e 



-i0 c 



Thus, 







(A.41) 



(A.42) 



gi(q,u n ) -gi(q,w n ) 

where gi(q, uj n ) is the Fourier transform anomalous correlator gf (x, r) = g\ (x, r). In 
what follows we denote the eigenvalues of G(q,u n ) as c±(q,u n ) = gi(q,u n ) ± gi(q,u„). 
Next we express Sqj, in terms of the eigenvectors of G(g, uS) and to subsequently 
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integrate out the transformed HS fields. In the following we use a matrix notation 
where the summations over R, R' as well as over q, u n and multiplication by h^ 1 are 
implicitly understood. Hence, 

[St, 5, Q\ Q]=\ [St (J- 1 + G) H + Q^GQ - 

Let us make the replacement G = UCUt, so that 

s™ a [x\x,p\p] = \ r 

where X = U f S (Jft = wtu), and P = U+Q (P f = QtU). We shall obtain the 
generating functional in terms of these new sources. We can now integrate over the 
auxiliary fields X\ X. Since they are related to Ej, H by a global unitary transformation 
the measure of the functional integral is not affected. Thus we are free to shift 



X ] ( J 1 + C) X + P f CP + P ] CX + X f CP 



(A.43) 



(A.44) 



X -f X - (J- 1 + C) 1 CP and Jft Jft _ ptQ (J- 1 + C) -1 , which yields 



P f ,P 



-iXtCX-ipt( c-C— r— C )P 



-ipt( C-C 



CP 



(A.45) 



where Af' is an unimportant constant. Although we have been careful to treat the 
matrices in the last expression as non-commuting, all of them are diagonal and the last 
expression can be simplified to: 



P f ,Pl = Z[0,0] exp 



ipt_ 

2 1 



C 



JC 



p 



A(Q,q,iu n ) 



Thus, the system has two modes, whose propagation is described by A 
When expressed in Fourier components this propagator reads: 

A+(Q, q,u n ) 

A_(Q,g,o; n ) 

A ± {Q,q,u) = — r 

1 + J(Q)c±(q,u n ) 

The excitations are poles of the propagator: 

1 + J (Q)c± (q, i^n —* w± (g,Q))l + J(Q)4(q,u ± (q,Q)) = 0. 



C(l 



(A.46) 
-JC)- 1 . 

(A.47) 
(A.48) 

(A.49) 



Before we solve these equations in detail, let us first demonstrate the existence of 
the Goldstone mode, namely that there is a solution of (|A.49J) at w = for |Q| 
and q — > 0. To this end, we use (jA.39IA.40|) . which are a consequence of gauge 
invariance, together with (|A.36|A.38J) . to arrive at c_(g = 0,u = 0) = gx(q — 0,u — 
0) - g x (q = 0,lu = 0) = S'^(\£ C \ 2 )/(M L(3). Moreover, by virtue of equation (|AT6jl . 
c-{q = 0,u = 0) = S'^ \ 2 )/(M L(3) = -EnJn = ~[J(Q = = -l/(J*c). This 
implies that u = is a solution of equation (|A.49J) for |Q|, q — > 0, as anticipated, and 
identifies the eigenvalue labeled the dispersion of the Goldstone mode. We 

shall use the condition that c_(g = 0,u = 0) = — [J(Q = 0)] _1 in the calculation of the 
full dispersion of this mode below. But before we proceed into this calculation we need 
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to know more about the correlation functions c±(k,u) for general k,u. Let us for a bit 
return to real space; we first notice that 

2 r -, 

c+(x,t) = gi(x,r) + gi(x,r) = -- \A B po (cos9(x,r) 0080(0, 0)) sG - ip 2 \ , (A.50) 

2 

c_(x,r) = gi (x,r) -gi(x,r) = ~^A B po (sin 9(x, r) sin 9(0, 0)) sG . (A.51) 

To make contact with the notation of Ref. [66J, we introduce $ = 9/(3, (3 2 = ff/K, 
x = (v s t,x), and [i = \/ A B PoCc/^ v s — \/A~bPo Jzc^c/^s- The excitations of such a 
sine-Gordon model are the solitons and anti-solitons. Furthermore, f3 2 < 4ir (i.e. K > 4) 
the model also has breather excitations |4TJ1 ITT] , which are soliton-anti-soliton bound 
states, and whose zero-momentum energies ("mass" gaps) are |4TJllT7j : 

A n = 2A, sin(7rpn/2), (A.52) 

with p = (3 2 / (8n — (3 2 ), n = 1, 2, ... , [p^ 1 ] — 1, and A s ~ jj, 2 - 1 /^ is energy gap of the 
soliton (see equation ([20)1 ). The operators sin/5$(x) and cos/?$(x) have zero conformal 
spin and therefore can only create breathers (if they exist) or pairs of solitons and 
anti-solitons. Furthermore, since under charge conjugation |4Uj . 

C- l $C = (A.53) 

it follows that 

C -1 (cos/3$)C = +cos/3$, (A.54) 

C~\an/3^)C = -sin/3$. (A.55) 

Thus the cosine operator can only have poles corresponding to breather states that are 
even under conjugation, whereas the sine can only have poles for breathers that are 
odd under conjugation. Taking all these facts into account we can write the spectral 
function of the correlators of these two operators |4T)] : 

p^(q,u)=2nZ £ w n 5(u 2 - v 2 s q 2 - A£) + pf^(q, lo), (A.56) 

n odd 

p™(q,u;) = 2nZ £ w n 5(u 2 - v 2 s q 2 - A^) + p^ s (q, u), (A.57) 

neven 

where p S ss L,cos (q,u)) describes the continuum of soliton and anti-soliton excitations. We 
have also used that the solitons with even quantum number n are even and those with 
odd n are odd under C [10]. To make further progress, we assume that the behavior 
of the above correlation functions is dominated by the lowest energy pole (single-mode 
approximation, SMA), i.e. the lowest breather Ai = 2A S sin(p7r/2). Hence, the retarded 
correlation function: 

z 

(u + ley — v z s q — (A 1 /h) 2 

where e — > + . Since c_(g = 0, uo = 0) = — [J(Q = 0)] 1 = —l/(Jzc), this fixes the 
residue z_ = Af / (JzcTi 2 )- Hence, the dispersion of the Goldstone mode can be obtained 
from equation (|A.49J) (we consider a square lattice below): 



Af J2j=y,z COS Qj a 

ui 2 _(qM)-v 2 q 2 -(A 2 /h 2 



0, (A.59) 
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so that u 2 _ {q, Q) = {v s q) 2 + ±A? E j=y , z (1 - cosQjO) /h 2 ~ (*; s g) 2 + (A 1 a/2^) 2 Q 2 , the 
last expression being valid for |Q| <C a -1 . However, we emphasize that note that as 
K — > +00 (i.e. the bosons become weakly interacting) the number of breather poles 
proliferate and the SMA may break down. Going beyond the SMA requires some 
knowledge of the spectral weights w n . Instead, we show the same form of the dispersion 
for the lowest-energy Goldstone mode can be also obtained from the self-consistent 
Harmonic approximation described in section 13.21 

Using similar methods we can obtain the energy dispersion of the longitudinal mode, 
which is the solution of the equation: 

l + J(Q)cJ(g,w + (g,Q)) = 0. (A.60) 

We again resort to the SMA, but taking into account that this time the pole is given 
by the lowest energy breather that is even under C, i.e. A2 = 2A S sin(7rp). We first fix 
the pole reside z + by studying the limiting behavior of c(q, lj) as q, uj — > 0. This can be 
obtained from (jA.36IA.38|) . along with (jA.39IA.40|) . which imply: 

c + (q = 0, u = 0) = 9l (q = 0, u = 0) + g^q = 0, u = 0) (A.61) 

1 r S>Mf)+2\tfS"(\Z c \ 2 )}=^^ (A.62) 



h d 2 £,n H ( 1 



(A.63) 



2(Jz c ) 2 d^ 2 c z c J \8K-l 
where we have used gauge invariance to set 6 C = so that |£ c | = £ c , together with 

d^ 2 8K-V K J 

which can be directly obtained from the expressions for E s q and ip c in section I3~T1 Hence, 
assuming that (SMA): 

c+(q,uj) ~ Z -±~ 7-r /M2 , (A.65) 

we find that z + /(A 2 /h) 2 = {zcJ)' 1 [1/ {8K - I)]. If, within the SMA, we solve 
equation (jA.60|) . we arrive at: 

u>l(q, Q) = A+ + v 2 q 2 + ^ (g^y) F(Q), (A.66) 

where A^ = A 2 , (8K — 2)/(8K — 1) is the energy gap of the longitudinal mode at 
(g, Q) = (0, 0), and F(Q) = a~ 2 £ t (l - e***) /2. 

At the transition temperature, T = T c , the order parameter ip c = £ c = and 
the matrix Jr^r/ + Gr_r/ becomes singular. For £ c = 0, the anomalous correlators 
gi + = g±~ = 0, and gf~ = + = g\, which leads to: 

^r-r'^( x — x')8(t — t') + gi(x — x' , t — t') = 0. (A. 67) 

Upon multiplying by Jr»_r, summing over R, and integrating over x and r, we obtain: 

/•ft/3 pL 

E J *-W dr / dx gi (x,T) = ZcJg?(q = 0,u = 0;T c ) = -l. (A.68) 
^ J0 Jo 
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This is precisely the same condition for the BEC temperature derived in section 13.1.21 
The expression for g^(q,u), which is needed to obtain T c explicitly is evaluated in the 
appendix that follows this one. 

Finally, to make connection with the discussion of section 13.1.31 we explain that 
the identification of the mode labeled as '— ' with the Goldstone mode and the one 
labeled as '+' as the longitudinal mode is indeed quite natural. If we perform a 
gauge transformation such that 6 C = 0, the matrix U of equation (jA.41|) produces 
the transformation: 

S&ix, t) = ±= [6£b.(x, r) - 6£ R (x, r)] , (A.69) 

Let us consider (small) fluctuations of the order parameter of the form considered in 
section EHlfl 6£{x, r) oc 5^ c (x, R, r) = * c (a;, R, T )-ip c = bPc + r] c {x, R, r)] e i6c -ip c ~ 
r) c (x, R, r) + i6 c (x, R, r), and S^(x, r) the complex conjugate. Hence 5^\x,t) oc 
i9 c (x, R,r) and 5^ (x, r) oc r] c (x, R, r). Thus we see that 5^ (x,r) is related to 
fluctuations of the phase of the BEC, whereas 5^ \x, r) is related to fluctuations of its 
amplitude. 

Appendix B. Finite-temperature phase susceptibility 

In this appendix we compute the Fourier transform of the retarded phase-susceptibility 
of a ID system of interacting bosons. This correlation function is to leading order the 
one-body boson Green's function: 

9i (x,t,T) = -~<T k R (a;,r)^(0)l) ~ ~ J p A B (e ie ^e w ^) . (B.l) 



n 



h 



In the above expression, we have replaced the boson field by its bosonized form to 
leading order: 

^(x,r)^pl /2 A]i 2 e w ^ T \ (B.2) 

where As = Ab{K) is a non-universal prefactor. To compute this correlation function, 
the (gaussian) action (|A.22|) must be used. Either by direct computation using 
functional integrals or by means of a conformal transformation of the correlation 
function at T = [37], one obtains the following result: 

f T , 1 , fK ~ ( (ttT/^pq) 2 \* 

9i{x,r,T) = --PoAb(K) - , B.3 

h \sm [irl [y s T + ix)/nv s \ sm [7ri (t> s r — tx)/hv s \J 

In order to obtain the retarded version of this correlation function, we shall use the 
relationship (see e.g. [17 \): 

g?(x, t, T) = -20{t) Im g[(x, t, T) . (B.4) 
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where gf(x,t,T) is the time-ordered (a la Feynman) propagator. The latter can be 
obtained from ()B.1|) by analytical continuation: r = it + e(t), where e(t) = sgn(t)e 
(e -> 0+). Hence, 



(TrT/hv s po) 



® l ^ ' ' ^ ^PoAb( ) ^ sinh. [ttT 1 (a; + v s t e /^,f s ] sinh[7rT(a; — v s t e )/hv s ] 
where t t = t — ie(t). We next exponentiate the above result and use that 

R(x, t) = sinh [tiT(x + v s t — iv s e(t))/hv s ] sinh [tiT^x — v s t + iv s e(t))/hv s ] 
= sinh [irT(x + v s t)/hv s ] sinh [irT(x — v s t)/hv s ] 



+ ie 



ttT 

T 



sinh 



/2nTt\ 



1 



0(e 2 



(B.5) 



(B.6) 



(B.7) 



which follows upon expanding in powers of e, as well as the identity ln(a + it) = 
In \a\ + m6{— a) (the branch cut of the logarithm is put on the negative real axis). Thus, 
the imaginary part of gf(x, t, T) does not vanish provided that ReR(x, t) < 0, that is, if 
sinh [nT(x + v s t)/hv s ] sinh [tiT(x — v s t)/hv s ] = — sinh [7rT£ + /hv s ] sinh [7rT£_ /hv s ] < 0, 
i.e. for > (£± = v s t ± x). Therefore, 



Img((Z + ,^,T) = -p A B (K) 



ttT 

hv s po 



D(K) 



' ttT ' 

hv s p 



e -^lnRei?( ?+ , ? _) Imexp 
1/4K 



17T 

AK 



(B.8) 



sinh [irT^+fTiVs] sinh [irT^_/hv s 

where D(K) = p Ab(K) sin (4^) Since we are interested in the Fourier transform 
of gi(x, t, T), it is convenient to express it in terms of the imaginary part of g[(x, t, T): 



g?{k,u,T) 



-2 
-2 
-2 



+00 



dt 



-00 
+00 



dt 



dt 



+00 

-00 
+00 

-00 

+v s t 

-Vst 



dxe 



dx e 



-ikx+iujt „R 



-ikx+iut 



g?(x,t,T) 
Im g F (x,t,T) 



dx e 



-ikx+iujt 



Im g((x,t,T) 



- f °>+ / + °V e^- +fc -^)/ 2 Im^f(e + ,^,T), 
v s Jo Jo 



(B.9) 
(B.10) 

(B.H) 
(B.12) 



where we have used that dtdx = d£, + d£- /2v s and have introduced k± = u/v s ± k. In the 
third step we employed that Im g[ (x, t, T) is only non-zero for — v s t < x < v s t. Finally, 
we have expressed the integral in terms of "light-cone" coordinates £-t to be able to 
benefit from the separability of Im g[ in terms of these coordinates. Hence, 



g?(.k,u,T) 



D(K) ( ttT 



where (see e.g. [T7] 



HVsPoJ 

e ik ^ 2 di 



f 



hv s ( uj , ... 



f 



T 



uj 

— q 

Vs 



2«? 



2T 



B 



8K 



.Kv s q 
1 — — , 1 
47rT' 



K 



AK 



(B.13) 



(B.14) 
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B(x, y) being the beta function, B(x, y) = T(x)T(y)/T(x + y). In the calculation of the 
critical temperature we need the value of gi(k, uu, T) for q = and u = 0: 

Appendix C. Self-consistent Harmonic Approximation calculation 

In this appendix, we fill in some of the details of the SCHA calculation. First, note that 
in equation we can rewrite 

(S-Sv)y = ~ G v (q, Q, UnjGv 1 (q, u n ) + (S- mt ) v - const.. (C.l) 

q,Q,d>n 

Then, the variational approximation to the free energy becomes 

F'[Gv} = -\ E lnG v (g,Q,^ n ) + ^ ]T G v (q,Q^n)G^(q,uj n ) + (S int ) v , (C.2) 

q,Q,U! n q,Q,U! n 

where 

(S- mt ) v = -A B p JL(3 E Re exp [G v (0, R - R', 0) - G v (0, 0, 0)] , (C.3) 

(RR') 

and 

G y (x, R - R, r) = — i— £ e iQ-(R-R')+^-i- n Gv(ffj q ^ n)) (C 4) 

In evaluating (5 , int ) v , we can take advantage of point group symmetry of the lattice in the 
perpendicular directions and replace G v (0, R — R', 0) — G v (0, 0, 0) inside the exponential 
by the more symmetric expression j- J2w=n+t G v (0, R — R', 0) — G v (0, 0, 0), where t are 
the unit lattice vectors in the perpendicular directions. Then, the outer S/rra just 
gives a factor of z c . Taking the variational derivative of F, we get: 



1 1 ZcAb ^ K)po JF(Q) e^^^'^'< m ' )GAq, ' Q '^ n) , 



(C.5) 

G Y (q, Q,u n ) G {q,u) n ) 2h 
with the free phonon propagator: G () 1 (q, u n ) = uj\ + v ^q 2 (v\\ = v s ), and 

F(Q) = J2t (l ~ eQ *) • This single self-consistent equation for G v (q, Q, u n ) can 
be solved by rewriting it into two equations: 



n = —(Poa) (^) V „ e ^OT^- F(Q)GvfeQ ^ } , (C.6) 



2K 



G;\q, Q, u; n ) = — k + (^||^) 2 + v\ F(Q)] . (C.7) 

Substituting equation (jC7|) into (|C.6|) . we get the equation for the variational parameter, 
the perpendicular phonon velocity v±: 

i- Hi _ 2q2 v F(n) ^ K (cr\ 
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The Matsubara sum can be done to give: 



7TV\\ 



In -. ^t~, — = -zrI(T, v ± ) (C.9) 



2K 



with 



r i V 2 

2a 2 _ F(Q) , + W 1 F (Q). 



We first study the variational equation (|C.9Jl at zero temperature. At T = 0, 
coth(x/2T) — > sign(x). To perform the summations over momentum, we take the 
continuum approximation for both the parallel (q) and perpendicular (Q) momenta, 
with a cut-off for the parallel momentum to be A = fiio/vu (fi\D is the ID chemical 
potential of an isolated tube). Thus, with J2q — > a 2 !Z{?/ a , and j^J2 q — > 

I— A we get: 



/( 0)W 0) w J- 
7ru|| 

where v± = v±(T = 0) and 



In 



k 4 j a J 2 



(C.ll) 



B = — r dQxdQ y F(Q) lna 2 F(Q) « 0.836477. (C.12) 

2 J-7T (271~) 2 

Putting equations (jC.llj) and (|C.9|) together then gives equation (|59*j) of Section 13.21 
We have assumed here v° ± /a <C Otherwise, our starting point using the Haldane 
harmonic fluid construction is not valid in principle. 

When K — > 00 (towards the non-interacting Bose gas limit), Ab{K) ~ 1, and using 
the Lieb and Liniger solution [T] we can show that fim ~ v\\p 7i/K. Thus, asymptotically 
for large K, we get: 

J!L~HiJ£ > SZ. (C.13) 

Thus, f°/a < yUi£i if and only if z c J < p, lD . But since for the nearly free Bose gas, 
\l oc 1/ i^ 2 is small, this demands a very hopping amplitude J. 

In equation ()59|) . there is a singularity as K — > 1/4, which is a signature that for 
K < 1/4 the Josephson coupling becomes irrelevant in the RG sense. As this regime is 
beyond the reach for bosons with a Dirac-delta interaction, we will not consider it any 
further. 

We now look at the variational equation (|C.9|) at finite temperature. First, defining 
v±{T c ) — v<l, we rewrite (|C.9|) as: 

/ \ 1-1/4K 

[jt ) = ^ [I(T e , vl) - 7(0, . (C. 14) 

Following Donohue[65j, we approximate crudely coth(x) « 1/x for |x| < 1, and 
coth(x) ~ 1 otherwise. Next, we develop a series expansion in v\jaT c < 1 for the 
integrand in I(T c ,v^_) — 1(0, Vj_). It will turn out that at the transition T = T c 



- CI 
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actually finite (i.e. first order transition) but less than Vj_, and the expansion can be 
justified a posteriori, at least for K not too close to 1. With the help of Mathematica 
to perform the expansion and to evaluate the resulting integrals, we find: 

V „ (J(T C , vl) - 7(0, «D) « + c - i In ^ + £ c 2n f^f* , (C.15) 

where c_ a « 0.677473, c « 0.133129 - 1/tt, c 2 = 5/96vr, c 4 = 21/2560vr, 
c 6 = 845/344604vr, c 8 = 9415/9437184vr, c 10 = 112203/230686720tt, c 12 = 
233695/8724152327T. We can show that the series converges quite rapidly for Vj_/ aT c < 1, 
but because this expression appears in an exponential, we need quite a few terms to 
get accurate trends, especially when K gets close to 1. Introducing this expansion into 

cm, 



v_l = (zcaTA 1/4K -&[e-i$£ -h»+£^i 

< " V < ) e 



(C.16) 



This equation is solved graphically for (T c , Vj_) (the right hand side of f)C. 16|) is an s- 
shaped curve that intersects the straight line of the lef hand side of (jC.16|) at 0, 1, or 2 
points, as a function of T c , and Vj_; T c is, by definition, when the intersection is at one 
point only. We use Mathematica's routine FindRoot to bracket this point to the desired 
accuracy.) 

We thus find the following asymptotic results for K — > oo: 



rri 

" ' - 0.610 + 0.69871", (C.17) 



% w 0.376 + 0.274/7T. (C.18) 

Vj_ 

Both asymptotes appear to hold good when K > 4. Dividing these equations, we get 
that aT c /vj_ ~ 1.86K + 0.270 + • • •, thus justifying our expansion in Vj_/aT c , at least for 
K > 4. Combining with the asymptotic results for 70 as K —>■ 00, we get the results 
quoted in section I3~2l 

Appendix D. Derivation of RG equations 

In this appendix we show how to compute the RG flow to second order in the Josephson 
and Mott potential couplings. To this purpose, it is convenient to introduce the chiral 
fields 0r+(x) and 0r-(x), which are implicitly defined as follows: 



k(z) = ^ [<h+(x) + R -(:r)] , (D.l) 
0k(x) = [<h.+(x) - <Pn-(x)] . (D.2) 



The effective Hamiltonian, equation thus takes the form: 
H e s = ~^-J2jdx (9 a .0R + (x)) 2 + (d x (j)R-(x)) 2 
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+ 9uYl I dx cosV^(0r+(x) + <pn-{x)) + g F E / dxd x (j) n+ (x)d x (j) n >. 

R " "' 

Vr+W - 0R'+(^) 



(R,R> 



c?x COS 



(R,R) 
^R-(z) 



>R'-K X 



2VK 



2VK 



(D.3) 



We have added here a term proportional to g F , which is not present initially, but which 
is generated at second order in gj after integrating out short-distance degrees of freedom 
(see below). It represents a density interaction between nearest neighbor ID systems. 



We next introduce the chiral vertex operators Vg(R,z) 



3 i/30 R+ (2) . and 



• e i/30R_(z) anc | cl e f| ne z = VgT _|_ ^ and 2 = t> s r — ix, and gPz = v s drdx = 
dzdz/2i. To explicitly display the scaling dimensions of the various operators in H e g, 
we introduce the dimensionless couplings: 

g u = Tralg u /v s , gj = 7ia 2 gj/v s g F = ng F /v s , gK = ng K /v s . (D.4) 

Thus, the interactions are described by 

r (Pz — r (P z — 

Smt = gi<Yl / TT" 90r + (^)50r_(^) +g F / IT- d(j) n+ (z)d(j) RI -(z) 



+ J^Yj^ ^(R, *) + (VF 

R, 



+ 



2-1/2A" 



a, 







(R,R) 



d 2 Z 



E / V" V-^(R,^)y+ 1 (R',0)l/_- 1 (R^)lCi (R'^) 

' Z7T — " 



2VK 



2\f~K 



2VK 



+ (Vk 



(D.5) 



The RG is performed in real space on the perturbative expansion of the partition 
function Z(a ) = Z (Texp [— S{ nt \). We first consider an infinitesimal change of the 
short distance cut-off a — > a' = (1 + S£)a , where < W < 1. Subsequently, we 
integrate out short-distance degrees of freedom in the range a < \z\ < a (l + 5£)a . 
Following Ref. |S7} . we use operator-product expansions (OPE) to evaluate expectation 
values of products of operators at two nearby points in space and time. For z — > w 
(respectively, z — ■> w) : 



d(p R _(z)d(j) R ^(w) = 

<90r+(^)<90r+(w) = 
V+(R,z)V^(R,w) 
Vp(R,z)VZ p (R,w) 
d<f> R _(z)V p -(R,w)-. 
3<f> K+ (z)Vp + (R,w)-. 



(z — w) 2 

1 



+- : 



>R 



(z — w) 2 

1 

(z — w) 13 

1 

(z — w) 13 
-i/3 



+ : 



<90R H 



z + w 



z + w 



l + ifi(z- z)d4>-R + 
1 + i(3(z — w)d(f)-R,- 



: +••• 

: +••• 

z + w 
2 

z + w 



+ ■ 

+ 



z — w 
z — w 



(D.6) 
(D.7) 
(D.8) 
(D.9) 
(D.10) 
(D.11) 
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Therefore, we begin with by considering the perturbative expansion of the partition 
function at the new scale a' n 



Z((l + 6£)a ,{ 9l (£ + 6£)}) 



ZW({l + 6£)a ,{ gi (£ + 6£)}) = - (S int ), 
z^((i + 6£)a 0l {g>(£ + S£)}) = ^(SL) 



= (1 + 5£)a , 

1 + Z^((l + S£)a ,{g t (£ + S£)}) 
+ Z^((l + 5£)a , { 9i (£ + 5£)}) + ---, (D.12) 

where 

(D.13) 
(D.14) 

The partition function must be left invariant by the RG transformation, i.e.Z((l + 
5£)ao, {gi(£ + 5£)}) = Z(a Q , { gi (£)}), provided that the couplings are properly 
transformed after integrating out the short-distance degrees of freedom in range defined 
by a < \z — w\ < (1 + 5£)ao. 

At the lowest order (tree level), the couplings change only because of the explicit 
factors of a (z.e.the scaling dimensions of the perturbing operators) in (jD.5|) . Thus, 
the only couplings that change are: 

9u(£) = 9u(£ + 5£)(1 + 5£)- {2 - K ^ ~ g u {£ + 51) (1 - (2 - K{£)5£)) , (D.15) 
and 

1 



qjW = gj(t + se)(i + siyv-v 2 ^ ~ gj (t + 51) 



2 - 



2K{£), 



SC. 



(D.16) 



At second order, there are ten terms in the expansion of Z^[(l + 5£)ao, {gi(£ + 5£)}}. As 
the manipulations are standard but rather long, we will just look at two of them in detail 
to illustrate the procedure. To see what kind of terms are generated from the product 
of two operators, we shall use "center-of-mass" coordinates: u = z — w (respectively, 
u = z — w) and v — (z + w)/2 (respectively, v = (z + w)/2), and we split the integrals 
of the operator products appearing at second order as follows: 
<Pv 
27 



d 2 u „ , u «,„ , u u. 
^0 1 (v + -,v + -)0 2 (v--,v--] 



(l+5£)a <\u\ 



2vr 



d 2 u 
~2^ 



d 2 u 
~2^ 



Oi(v + -,v + -)0 2 (v--,v-- 



(D.17) 

L ao<| , u| ao<\u\<(l+5l)ao 

The terms generated at second order by the RG transformation stem form the second 
integral, where the relative coordinate is restricted to the infinitesimal range < \u\ < 
a (l + 5£). To O(gpgj), there are two possible operator pairings because the hopping 
term couples operators from neighboring chains: 



0(g F 9u) 



1 9f{£)9j{£) 
2! al 1/2K 

V + i (T, v - 

2\fK 



(R,R'),(T,T'> 



d 2 U f/K , / U,., , W, 

l(d(f> n+ (v + -)d ( f) W _(v + - 



a <\u\<(l+5e)a 



:D.18) 
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Focusing on the terms where R = T and R' = T', or R = T' and R' = T and using 
the OPE's (jD.10ID.ll|) with (3 = ±1/2 \/K, we arrive at the following expression: 



0{g F g u ) 



2 x 2 x 



1 9f(£)9j(£) 
2! a 2 - 1/2K 



AK(£) 



E 

(R,R') 



~2^r~ 



V^(R,v)V+^(R',v) 

2V7f 2\f~K 



x Fx (R,u)H (R» + (VK 

IsTE 2-Jk ^ 

g F (£)gj(£)5i rd 2 v~ 

(R,R'> 



(l+«)ao d U 
a U 



(D.19) 



4K(£)ao _1/A ' 



2tt 



V\_ (R, (R, v)VZ_l_ (R, (R', v) 

2%/~K 2VK 2VK 2VK 



(D.20) 



where the first factor of two stems from the two terms of O(gFgj) and the second from 
the two possible pairings of R, R' and T, T'. Thus we see that to second order in the 
couplings the coupling gj renormalizes as follows: 



gj(£)=gj(£ + 5£)\l 



(2 



-)gj{t) + 



g F (i)gj(i) 



In the limit 5£ —>■ 0, this leads to the following differential equation: 



81. 



dg.j 
d£ 



2K ]9J + 



9j9f 
AK ' 



(D.21) 



(D.22) 



Other terms are dealt with in a similar fashion. It is also worth noticing that terms of 
0{g 2 Kl gFgx, gxgj, 9k9j) need not be taken into account because at each RG step we 
set gx = by properly renormalizing K. Thus gx{£ + 8£) = 0, but new terms of the 
form of the operator proportional to g% are generated at 0(g 2 ,g 2 ). In order to get rid 
of them, we perform an infinitesimal canonical transformation. To understand this, let 
us consider a simplified version of the above model, where the chain index R is dropped, 
and the Hamiltonian reads: 



H = 
+H. 



Air 



dx 



(d x ci> + (x)f + (d x <p.(x)f 



+ 



hv s 5g 
2ir 



K 



+ 



2VK 



d x <j) + (x)d x (f>-(x) 

(D.23) 



We assume the dimensionless coupling 8g -C 1. The precise functional forms of H± are of 
no importance in what follows. Next consider the following infinitesimal transformation: 

0+ = <P+ ~ y^-, 



(f>- = <P~ ~ y</>+- 

To 0(6 g) one can show that this is a canonical transformation which, to 0(5g 2 ) 
the first term of the Hamiltonian to the diagonal form: 



Air 



dx 



(d x (f) + (x)f + (d x <p-(x) 



+ 0(8g 2 ). 



(D.24) 

(D.25) 
brings 

(D.26) 
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However, the combinations: 



+ 0- = (l-y)(0 + + 0-), (D.27) 
_ 0_ = (1 + & (} + - 0_) . (D.28) 



2 

Since these combinations enter H + and H_, it means that K must be renormalized as 
follows: 

K=U-^j K~(l-6g)K (D.29) 

Returning to the problem of interest, for the array of coupled ID Bose gases we have 
that 8g = g K (£) = (gj/K(£) - g 2 v {l)K{l)) 51 Hence, 

K{£) = (1 - 5g)K{i + 5t) = [l - (ng 2 j/K(£) - g 2 v {£)K{£)) 5i\ K{£ + 5£), (D.30) 
which can be turned into the following differential equation: 

d -^ = ng)- g 2 v K 2 (D.31) 
With these considerations, the equations given in section 14.11 can be obtained. 

Appendix E. Spin chain mapping 

In this Appendix, we outline the mapping from the continuum double sine-Gordon 
model to a spin chain model. We take the continuum limit of the spin chain, and use 
the standard spin-operator bosonization formulas [T71 l4"U]: 

— S z m -> S z (x = ma ) = -d x (f> + — (-1)*/°° cos20(x), (E.l) 
ao vr 7rao 

-S+ - S + (x = ma,) = K -^=e- ie ^\ (E.2) 

where ao is the lattice spacing of the spin chain. Using this mapping, the continuum 
double sine-Gordon mean field model (|78j) becomes the spin chain model of (fT^j) in 
Section 14.21 The mean field self-consistency condition is still as before (|18jl. For 
simplicity, we take the isotropic Heisenberg chain. Then v s = v\\ = nJ a /2, and 
K = 1/2. Now K = 1/2 is not a physically valid regime for bosons interacting with a 
contact potential. However, from the RG approach above, there are no indications of 
a phase transition from K = 1/2 to the physical K > 1 regime, and we expect that 
at least qualitatively, the results derived here should be valid for K > 1. Moreover, 
strictly speaking, there is a umklapp operator ~ cos40(x) that is marginally irrelevant 
at K = 1/2, but since it is less relevant than the Mott and Josephson terms we drop it. 

Appendix F. Derivation of quantum phase model 



Consider a 2D array of finite tubes described by the Euclidean action shown in 
equations (jjDj) and (|5Uj) . We describe the tubes as ID systems with open boundary 



Quasi-ID interacting Bose gases 



53 



conditions |3T?1 13*7] . In such a finite systems, the phase field operator consists two 
terms: 



Or(x,t) = 9 0R (t) + R (a;,r) 



(F.l) 



where 6r(x, r) describes phase fluctuations (phonons) of wave number g ^ 0. The 
operator # r is canonically conjugate to the number operator iV R , i.e.: [0 O r, N-r/] = 
*<5r,r'- Introducing (|F.1J) into pHjl. the quadratic part of S reduces to 



2nv s ^ Jo 



(d T R (x, r)) 2 + v s (d x e R (x, t))' 



(F.2) 



The hopping term ([50)1 . however, couples ^or( t ) and ©r(^,t) in a non-linear fashion: 



Sj [#or, @r] = T^^bPo ^2 I dr 

(R,R'> 



da; 



e i(e O R(r)-0 OR /(r)) e i(e R (*,r)-e Il ,(x,r)) + c c 



(F.3) 



In order to obtain an effective action in terms of ^or(t) only we need to integrate out 
6r(x, r). We shall do this perturbatively assuming that J/fiio is small. Thus, let us 
define: 



-Sf [0o] = ( g -S J ) e 



e = exp 



■(Sj) e + ~({Sj-(Sj) e )*)e + 



(FA) 



where we have used the cumulant expansion and employed the following notation: 
_ J [d&] A[Q] e^gj 

" / [ d 0] e-*ie] • (F ' 5) 

for a given functional of G, A[G]. The lowest order contribution to Sj [9q] comes from 
the term (Sj)q, which yields: 



s e r ~ 



-AbPoL £ e"< e » rVcos[0 O R(T)-0oR'(T)] (F 

(R,R'> 



.6) 



Using that (O R (0))e = 1 ln(L/a ), where a is the short- distance cut-off, the above 
expression leads to: 



S 



eff 



1a b ( Po L)(-) " J2 rdrco8[e OIl (r)-eow(r)] (F.7) 

" Vao/ (R,R'} 



(aoPo 



,1/2X 



Assuming that each ID tube is described by the Lieb-Liniger model, Ab 
(K/-k) 1 I 2K |37j, where a ~ %v s / ' [i (i.e. approximately the healing length, £). Hence, the 
effective hopping energy of the quantum-phase model is Ej = Ej(N ) ~ J(N ) 1 " 1 ^ 2K , 
where we have used that p L = N . The complete effective action reads: 



S eS [9 ] 



hi3 



2E c Jo 

h 



dr (d T 9on(r)) 

rh(3 



rtip 

E / dr cos[fl 0R (r)-£ 0R ,(r)]+O(J 2 ), (F.8) 



<R,R'> 



Quasi-ID interacting Bose gases 



54 



where the charging energy Eq = hirv s /KL. This euclidean action corresponds to the 
following Hamiltonian (we restore the chemical potential term): 

R R 

-Ej J2 cos[6 0R (t)-6 0R ,(t)]+O(J 2 ). (F.9) 

(R,R'> 

This is the quantum-phase model used in section 15.11 
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